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Abstract 

We study the evolution of a particle system whose genealogy is given by a supercritical contin- 
uous time Galton- Watson tree. The particles move independently according to a Markov process 
and when a branching event occurs, the offspring locations depend on the position of the mother 
and the number of offspring. We prove a law of large numbers for the empirical measure of indi- 
viduals alive at time t. This relies on a probabilistic interpretation of its intensity by mean of an 
auxiliary process. The latter has the same generator as the Markov process along the branches 
plus additional jumps, associated with branching events of accelerated rate and biased distribution. 
This comes from the fact that choosing an individual uniformly at time t favors lineages with more 
branching events and larger offspring number. The central limit theorem is considered on a special 
case. Several examples are developed, including applications to splitting diffusions, cellular aging, 
branching Levy processes. 

Key words. Branching Markov process. Branching diffusion. Limit theorems. Many-to-One formula. 
Size biased reproduction distribution. Size biased reproduction rate. Ancestral lineage. Splitted 
diffusion. 
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1 Introduction and main results 

We consider a continuous time Galton- Watson tree T, i.e. a tree where each branch lives during an 
independent exponential time of mean 1/r, then splits into a random number of new branches given by 
an independent random variable (r.v.) v of law (p^, k £ N), where N = {0, 1, . . .}. We are interested 
in the following process indexed by this tree. Along the edges of the tree, the process evolves as a 
cadlag strong Markov process (Xt)t>o with values in a Polish space E and with infinitesimal generator 
L of domain D(L). The branching event is nonlocal: the states of the offspring are described by a 
function (F^ k \x, O), . . . , F^ k \x, 0)), which clearly depends on the state x of the mother just before 
the branching event, and of the number v = k of offspring ; besides, the randomness of these states 
are modelled via the random variable G, which is uniform on [0,1]. Finally, the new born branches 
evolve independently from each other. 

This process is a branching Markov process, for which there has been a vast literature. We refer 
to Asmussen and Hering [2], Dawson [T2] and Dawson et al. [IB] for nonlocal branching processes 
similar to those considered here. Whereas the literature often deals with limit theorems that consider 
superprocesses limits corresponding to high densities of small and rapidly branching particles (see e.g. 
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Dawson |15j , Dynkin [TO] , Evans and Steinsaltz [23J ) , we stick here with the discrete tree in continuous 
time which we aim characterizing. 

Let us also mention some results in the discrete time case. Markov chains indexed by a binary 
tree have been studied in the symmetric independent case (see e.g. Athreya and Kang [4J, Benjamini 

(2) (2) 

and Peres [8]), where for every x, F^'(x,e) and Ff>(x,0) are i.i.d. A motivation for considering 
asymmetric branching comes from models for cell division. For instance, the binary tree can be used 
to describe a dividing cell genealogy in discrete time. The Markov chain indexed by this binary tree 
then indicates the evolution of some characteristic of the cell, such as its growth rate, its quantity of 
proteins or parasites... and depends on division events. Experiments (Stewart et al. [S3]) indicate 
that the transmission of this characteristic in the two daughter cells may be asymmetric. See Bercu et 
al. [U] or Guyon [32] for asymmetric models for cellular aging and Bansaye [B] for parasite infection. 
In Delmas and Marsalle [17] a generalization of these models where there might be 0, 1 or 2 daughters 
is studied. Indeed under stress conditions, cells may divide less or even die. The branching Markov 
chain, which in their case represents the cell's growth rate, is then restarted for each daughter cell at 
a value that depends on the mother's growth rate and on the total number of daughters. 

We investigate the continuous time case and allow both asymmetry and random number of off- 
spring. Let us give two simple examples of this model for parasite infection. In the first case, the 
cell divides in two daughter cells after an exponential time and a random fraction of parasites goes 
in one of the daughter cell, whereas the rest goes in the second one. In the second case, the cell 
divides in k daughter cells and the process X is equally shared between each of the k daughters: 
Vj G {1, . . . , k}, Fj(x, O) = x/k. Notice that another similar model has been investigated in Evans 
and Steinsaltz [23J where the evolution of damages in a population of dividing cells is studied, but 
with a superprocess point of view. The authors assume that the cell's death rate depends on the 
damage of the cell, which evolves as a diffusion between two fissions. When a division occurs, there 
is an unbalanced transmission of damaged material that leads to the consideration of nonlocal births. 
Further examples are developed in Section [SJ 

Our main purpose is to characterize the empirical distribution of this process. More precisely, if we 
denote by N t the size of the living population V t at time t, and if (Xf) u& y t denotes the values of the 
Markov process for the different individuals of Vt, we will focus on the following probability measure 
which describes the state of the population 



This is linked to the value of the process of an individual chosen uniformly at time t, say U(t), as we 
can see from this simple identity: 



We show that the distribution of the path leading from the ancestor to a uniformly chosen individ- 
ual can be approximated by means of an auxiliary Markov process Y with infinitesimal generator 
characterized by: V/ G D(L), 



1 {N t >0} 

Nt 



uGVt 



J2Sx?(dx), t G R+. 



E 



' 1 {A f i>0} 

L N t 



u£Vt 



f(X?) =E[l {Nt>Q} f(X^)]. 




(1.1) 
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where we recall that r denotes the particle branching rate and where we introduce m = Ylk=l ^Pk 
the mean offspring number. In this paper, we will be interested in the super-critical case m > 1, 
even if some remarks are made for the critical and sub-critical cases. The auxiliary process has the 
same generator as the Markov process running along the branches, plus jumps due to the branching. 
However, we can observe a bias phenomenon: the resulting jump rate rm is equal to the original 
rate r times the mean offspring number m and the resulting offspring distribution is the size-biased 
distribution (kpk/m,k £ N). For m > 1 for instance, this is heuristically explained by the fact that 
when one chooses an individual uniformly in the population at time t, an individual belonging to 
a lineage with more generations or with prolific ancestors is more likely to be chosen. Such biased 
phenomena have already been observed in the field of branching processes (see e.g. Chauvin et al. 
[14] . Hardy and Harris [33] . Harris and Roberts [34] ). Here, we allow nonlocal births, prove pathwise 
results and establish laws of large numbers when Y is ergodic. Our approach is entirely based on a 
probabilistic interpretation, via the auxiliary process. 

In case Y is ergodic, we prove the laws of large numbers stated in Theorem 11.11 and 11.31 where W 
stands for the renormalized asymptotic size of the number of individuals at time t (e.g. Theorems 1 
and 2 p. Ill of Athreya and Ney [S]): 

W := lim N t /E[N t ] a.s. and {W > 0} = M > 0, N t > 0} a.s. 

t— >+oo 

Theorem 1.1. If the auxiliary process Y is ergodic with invariant measure tt, we have for any real 
continuous bounded function f on E: 

l im J2 f(X?) = 1 {W>0} f /(x)Tr(dx) in probability. (1.2) 

This result in particular implies that for such function /, 

lim E[/(Xf (t) ) | N t > 0] = / f(x)7r(dx), (1.3) 

where U(t) stands for a particle taken at random in the set Vt of living particles at time t. 

Theorem II .11 is a consequence of Theorem l4.2l (which gives similar results under weaker hypotheses) 
and of Remark 14.11 The convergence is proved using I? techniques. 

Theorem 1 1 . 1 1 also provides a limit theorem for the empirical distribution of the tree indexed Markov 
process. 

Corollary 1.2. Under the assumption of Theorem 

lim i^ t>0 } 5x u (dx) = 1 /iy>o> n{dx) in probability (1-4) 

1 uGVt 

where the space Mf{E) of finite measures on E is embedded with the weak convergence topology. 

We also give in Propositions 16.11 and 16.41 a result on the associated fluctuations. Notice that 
contrarily to the discrete case treated in [T7], the fluctuation process is a Gaussian process with a 
finite variational part. 

In addition, we generalize the result of Theorem II .11 to ancestral paths of particles (Theorem 1 1.3|) : 
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Theorem 1.3. Suppose that Y is ergodic with invariant measure ir and that for any bounded measur- 
able function f, lim^+oo Ea;[/(Yj)] = f E f(x)ir(dx), then for any real bounded measurable function ip 
on the Skorohod space O([0, T],E), we have 

lim {N *> 0} V ip(X^, t-T <s <t) = EJip{Y s , s < T)] ls wm in probability, 

1 u&Vt 

where, for simplicity, X" stands for the value of the tree indexed Markov process at time s for the 
ancestor of u living at this time. 

Biases that are typical to all renewal problems have been known for long time in the literature 
(see e.g. Feller [24J, Vol. 2 Chap. 1). Size biased trees are linked with the consideration of Palm 
measures, themselves related to the problem of building a population around the path of an individual 
picked uniformly at random from the population alive at a certain time t. In Chauvin et al. [14] and 
in Hardy and Harris [33], a spinal decomposition is obtained for continuous time branching processes. 
Their result states that along the chosen line of descent, which constitutes a bridge between the initial 
condition and the position of the particle chosen at time t, the birth times of the new branches form 
a homogeneous Poisson Point Process of intensity rm while the reproduction law that is seen along 
the branches is given by (kpk/m,k E N). Other references for Palm measures, spinal decomposition 
and size-biased Galton- Watson can be found in discrete time in Kallenberg [39], Liemant et al. [41], 
and for the continuous time we mention Gorostiza et al. |31j . Geiger and Kauffmann [29], Geiger [2HJ 
or Olofsson [49] . Notice that biases for an individual chosen uniformly in a continuous time tree had 
previously been observed by Samuels [52] and Biggins [T2J. In the same vein, we refer to Nerman and 
Jagers |47] for consideration of the pedigree of an individual chosen randomly at time t and to Lyons 
et al. [32], Geiger [27] for spinal decomposition for size biased Galton- Watson processes in the discrete 
time case. 

Other motivating topics for this kind of results come from branching random walks (see e.g. Biggins 
[33] , Rouault [5T] ) and homogeneous fragmentation (see Bertoin [TQ1 E] ) • We refer to the examples 
in Section 5 for more details. 

The law of large numbers that we obtain belongs to the family of law of large numbers (LLN) 
for branching processes and superprocesses. We mention Benjamini and Peres [5], and Delmas and 
Marsalle [T7] in discrete time, with spatial motion for the second reference. In continuous time, LLNs 
are obtained by Georgii and Baake [30] for multitype branching processes. Finally, in the more different 
setting of superprocesses (obtained by renormalization in large population and where individuals are 
lost), Englander and Turaev [20] . Englander and Winter [21] and Evans and Steinsaltz [23] have 
proved similar results. Here, we work in continuous time, discrete population, with spatial motion, 
and nonlocal branching. This framework allows to trace individuals, which may be interesting for 
statistical applications. Our results are obtained by means of the auxiliary process Y, while other 
approaches involve spectral techniques, and changes of measures via martingales. 

In Section 2, we define our Markov process indexed by a continuous time Galton- Watson tree. We 
start with the description of the tree and then provide a measure-valued description of the process 
of interest. In Section 3, we build an auxiliary process Y and prove that its law is deeply related 
to the distribution of the lineage of an individual drawn uniformly in the population. In Section 
4, we establish the laws of large numbers mentioned in Theorem 11.11 and 11.31 Several examples are 
then investigated in Section 5: splitting diffusions indexed by a Yule tree, a model for cellular aging 
generalizing [T7J and an application to nonlocal branching random walks. Finally, a central limit 
theorem is considered for the case of splitting diffusions in Section 6. 
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2 Tree indexed Markov processes 



We first give a description of the continuous time Galton- Watson trees and preliminary estimates in 
Section [2.11 Section [2.21 is devoted to the definition of tree indexed Markov processes. 

2.1 Galton- Watson trees in continuous time 

In a first step, we recall some definitions about discrete trees. In a second step, we introduce continuous 
time and finally, in a third step, we give the definition of the Galton- Watson tree in continuous time. 
For all this section, we refer mainly to |18[ [351 HO] ■ 

Discrete trees. Let 

+oo 

U = |J (N*) m , (2.1) 

m=0 

where N* = {1, 2, . . .} with the convention (N*)° = {0}. For u 6 (N*) m , we define \u\ = m the gener- 
ation of u. If u = (tii, • • • i u n ) and v = (v\, . . . , v p ) belong to U, we write uv = (m, . . . , u n , v±, . . . , v p ) 
for the concatenation of u and v. We identify both 0u and u0 with u. We also introduce the following 
order relation: u -< v if there exists w € U such that f = uro; if furthermore w ^ 0, we write u < v. 
Finally, for u and « in W we define their most recent common ancestor (MRCA), denoted by u A v, as 
the element w E IA of highest generation such that w <u and w ^ v. 

Definition 2.1. A rooted ordered tree T is a subset ofU such that: 

(i) e r, 

(ii) if v £ T then u < v implies u £ T, 

(Hi) for every u £ T , there exists a number u u £ N such that if v u = then v y u implies v ^ T, 
otherwise uj £ T if and only if 1 < j < v u . 

Notice that a rooted ordered tree T is completely defined by the sequence (i> u , u £ U), which gives the 
number of children for every individual. To obtain a continuous time tree, we simply add the sequence 
of lifetimes. 

Continuous time discrete trees. For a sequence (l u ,u & U) of nonnegative reals, let us define: 
Vu £ U, a(u) = l v and f3(u) = l v = a(u) + l u , (2-2) 

with the convention q(0) = 0. The variable l u stands for the lifetime of individual u while a(u) and 
(5(u) are its birth and death times. Let 

U = ^x[0,+oo). (2.3) 
Definition 2.2. A continuous time rooted discrete tree (CT) is a subset T o/U such that: 
(i) (0,0) 6 1 

(ii) The projection ofT onlA, T , is a discrete rooted ordered tree, 
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(iii) There exists a sequence of nonnegative reals (l u ,u G IX) such that for u £ T , (u,s) G T if and 
only if a(u) < s < j3(u), where a(u) and j3(u) are defined by ( QS| ). 

Let T be a CT. The set of individuals of T living at time t is denoted by Vt 

V t = {u G U : (u,t) G T} = {u G T : a(u) < t < 0(u)}. (2.4) 

The number of individuals alive at time t is Nt = Card(Vj). We denote by Dt the number of individuals 
which have died before time t: 

D t = Card{u G T : (3(u) < t}. (2.5) 
For (u, s) G T and t < s, we introduce u(t), the ancestor of u living at time t: 

u(t) = v if (v^u and (v,t) eT). (2.6) 

Eventually, for (u, s) G T, we define the shift of T at (u, s) by 9i u S \T = {(v,t) G U : (uv, s + t) G T}. 



Note that 0( UjS )T is still a CT. 



Continuous time Galton- Watson trees. Henceforth, let (f2, J-, IP) be the probability space on 
which we work. 

Definition 2.3. We say that a random CT on (f2, J 7 , P) is a continuous time Galton-Watson tree with 
offspring distribution p = (p^ ,fcsN) and exponential lifetime with mean 1/r if: 

(i) The sequence of the number of offspring, (^ u ,ii G IX), is a sequence of independent random 
variables with common distribution p. 

(ii) The sequence of lifetimes (l u ,u G IX) is a sequence of independent exponential random variables 
with mean 1/r. 

(iii) The sequences [y u ^ G IX) and (l u ,u G IX) are independent. 

We suppose that the offspring distribution p has finite second moment. We call 

m = ^kp k and q 2 = ^(k - m) 2 p k , (2.7) 

fc>o fc>o 

its expectation and variance. The offspring distribution is critical (resp. supercritical, resp. subcrit- 
ical) if m = 1 (resp. m > 1, resp. m < 1). In this work, we mainly deal with the supercritical 
case. 

We end Section [2. II with some estimates on Nt and Dt. To begin with, the following Lemma gives 
an equivalent for Nt- 

Lemma 2.4. For t G M+ ; we have 

E[N t ]=e rim - 1)t , (2.8) 

r 2l _/ e^" 1 " 1 )* +(? 2 (m- l)" 1 + m )(e 2r ( m - 1 )*-e r ( m - 1 )') if m ^ I . . 

E[iV * J "\ l + frt ifm = l. [ J) 

If m > 1 there exists a nonnegative random variable W such that 

lim — r-^-r = W a.s and in L 2 , (2.10) 

*-H-ooE[iVi] 

> 0} = {Vt > 0, N t > 0} a.s. and P(W > 0) > 0. 
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Figure 1: Continuous time Galton- Watson tree. 



Proof. The process {Nt, t > 0) is a continuous time Markov branching process, so that the expectation 
and the variance of Nt are well-known, see [5] Chapter III Section 4. Almost sure convergence towards 
W is stated again in [S] Theorems 1 and 2 Chapter III Section 7. Finally, since the martingale 
(Nte- r{m - 1]t ,t > 0) is bounded in L 2 , we obtain the L convergence (e.g. Theorem 1.42 p. 11 of 
[371). ■ 

We also give the asymptotic behaviour of Dt, the number of deaths before t. 

Lemma 2.5. If m > 1, we have the following convergence a.s. and in L 2 : 

D t 



lim 

t^+oo E[D t ] 



W. 



with 



and W defined by < TOfl|) . 



E[A] = [m-l)-\e 



1 r r(m—l)t 



(2.11) 



(2.12) 



Proof. First remark that (Dt, t > 0) is a counting process with compensator (Jq riV s <is, t > 0). We set 
AN t = N t - iV 4 _ so that dN t = AN t dD t . To prove ([2TTT]) . it is sufficient to prove that e- 7 "^- 1 )* I 4 
goes to a.s. and in L 2 , where It = (m — 1)-Dt — Nt- Since I = (It,t > 0) satisfies the following 
stochastic equation driven by (D t , t > 0) 



dZ t = (m-l-AJV t )dDt, 
we get that I is an I? martingale. We deduce that d(I)t = s 2 rN t dt and: 



(2.13) 



E[I 2 1 =l + E[(/)t] =1+?V /* e »-(m-l) ads = 1 + _i!_ ( e r(m-l)t_^ / 2 M) 

1 ' jo m — 1 

which implies the L convergence 

of e -r(m-i)t It to o. Besides, (e-^™" 1 ) 4 l u t>0) is a supermartingale 
bounded in I? and hence the convergence also holds almost surely. ■ 

Example 1. Yule tree. The so-called Yule tree is a continuous time Galton- Watson tree with a 
deterministic offspring distribution: each individual of the population gives birth to 2 individuals that 
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is P2 = 1 (i-e. p = S2, the Dirac mass at 2). The Yule tree is thus a binary tree, whose edges have 
independent exponential lengths with mean 1/r. In that case, W is exponential with mean 1 (see e.g. 
[5] p. 112). We deduce from Lemma YIM that, for t E R+, 

E[iV t ] = e r * and E[N?} = 2 e 2ri - e rt . (2.15) 

Notice that (|2.15|) is also a consequence of the well-known fact that Nf is geometric with parameter 
e~ rt (see e.g. §5] p.105). 

2.2 Markov process indexed by the continuous time Galton- Watson tree 

In this section, we define the Markov process Xj = (Xf, (u,t) E T) indexed by the continuous time 
Galton- Watson tree T and with initial condition p. Branching Markov processes have already been 
the object of an abundant literature (e.g. [21 El El EZ1 ES])- The process that we consider jumps at 
branching times (nonlocal branching property) but these jumps may be dependent. 

Let (E,S) be a Polish space. We denote by V(E) the set of probability measures on (E,£). 

Definition 2.6. Let X = (Xt,t > 0) be a cddldg E-valued strong Markov process. LetF= (F^ k \l < 
j < k,k E N*) be a family of measurable functions from E x [0, 1] to E. The continuous time branching 
Markov (CBM) process Xj = (Xf, (u,t) E T) with offspring distribution p, exponential lifetimes with 
mean 1/r, offspring position F, underlying motion X and starting distribution fj, E V(E), is defined 
recursively as follows: 

(i) T is a continuous time Galton- Watson tree with offspring distribution p and exponential lifetimes 
with mean 1/r. 

(ii) Conditionally on T, X % = (Xf,t E [0,^(0))) is distributed as (X t ,t E [O,/3(0))) with X dis- 
tributed as \i. 

(in) Conditionally on T and X® , the initial positions of the first generation offspring (X^ u y 1 < u < 
Vfjj) are given by (F^ (X^,^_ , 0), 1 < u < v§) where is a uniform random variable on [0, 1]. 

(iv) Conditionally on X® , v§, (3$ and (X™, u \,l < u < v§), the tree-indexed Markov processes 
( w > ^ ^(u,q(u))1') f or 1 — u — u % are independent and respectively distributed as 
Xj with starting distribution the Dirac mass at X™, u \ . 

For x E E, we define F X (A) = P(^4|Xq = x) for all A E F, and denote by E x the corresponding 
expectation. For fi E V(E) we set in a classical manner ¥^(A) = f E ¥ x (A)p(dx), and write E^ for the 
expectation with respect to P^. 

For u E T, we extend the definition of Xf when t E [0,a(u)) as follows: X™ = X^\ where u(t), 
defined by (|2.6p . denotes the ancestor of u living at time t. 

Notice that for u E T, (Xf,t E [0, f3(u))) does not encode the information about the genealogy of u. 
We remedy this by introducing the following process (A",i > 0) for u £U: 

v^u{t) 

This process provides the birth times of the ancestors of u, as well as their offspring numbers. Notice 
that it is well defined for all u E IA contrarily to Xf . Indeed, the state of u at its birth time, X^Ou)' * s 
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Figure 2: Continuous time Markov process indexed by the Galton-Watson tree of Figured 



well denned only for u £ T, since it depends on the state of the parent and the number of its offspring. 
For u £W, the process (A", t £ [0, f3(u))) is a compound Poisson process with rate r for the underlying 
Poisson process (Sf,t > 0) and increments distributed as log(v) with v distributed as p, stopped at 
its (\u\ + l)-th jump. 

In the sequel, we denote by X" = (X™, A") the couple containing the information on the position and 
genealogy of the particle u. 

2.3 Measure-valued description 

Let Bb(E, R) be the set of real-valued measurable bounded functions on E and Mf{E) the set of finite 
measures on E embedded with the topology of weak convergence. For p £ M.f(E) and / S Bi,(E,~$l) 
we write (p, f) = J E f{x)p(dx). 

We introduce the following measures to represent the population at time t: 



where Vt has been defined in (|2.4|) . Note that (Z t , f) = J2 u eV t fi^-t)- Since X is cadlag, we get that 
the process Z = (Z t ,t > 0) is a cadlag measure- valued Markov process of D(R+, M.p(E)). 

Following the work of Fournier and Meleard [26J , we can describe the evolution of Z in terms of 
stochastic differential equations (SDE). Let p(ds,du,dk,d0) be a Poisson point measure of intensity 
r ds ® n(du) 8) p(dk) (g> d9 where ds and dO are Lebesgue measures on R + and [0,1] respectively, 
n(du) is the counting measure on U and p(dk) = J2ie?$ Pi$i(,dk) is the offspring distribution. This 
measure p gives the information on the branching events. Let L be the infinitesimal generator of X. If 
C fe ' u (lR + x E,R) denotes the space of continuous bounded functions that are C 1 in time with bounded 
derivatives, then for test functions / : (t, x) h4 ft(x) in C^'°(IR + x E, M) such that Vt G R+, ft € D(L), 



Z t= 5 (u,x?), and Z t = ^ 



U 



(2.16) 




we have 



(Z t Jt) =fo(X$) + 



L 




huevs-} E fs(Ff\XL ,0))- / S (X«_ ) d«, dfe, d9) 





(2.17) 



where W/ is a martingale. Explicit expressions of this martingale and of the infinitesimal generator 
of (Zt,t > 0) can be obtained when the form of the generator L is given. 

Example 2. Splitted diffusions. The case when the Markov process X is a real diffusion (E = R) 
is an interesting example. Let L be given by: 

Lf(x) = b(x)f'(x) + ^f(x), (2.18) 

where we assume that b and a are bounded and Lipschitz continuous. In this case, we can consider 
the following class of cylindrical functions from into M. defined by 4>f(Z) = <f)((Z,f)) for 

/ G C 2 (R,M) and <f> £ C%(R) which is known to be convergence determining on 'P(yVip'(lR)) (e.g. [T5] 
Theorem 3.2.6). We can define the infinitesimal generator C of (Zt)t>o for these functions: 

C4f(Z) = d<j> f (Z) + C 2 <l>f{Z), (2.19) 

where C\ and £2 correspond to the branching and motion parts. Such decompositions were already 
used in Dawson [T5] (Section 2.10) and in Roelly and Rouault [50] for instance. The generator C\ is 
defined by: 

C^f{Z)=r f L((Z,/) + E/(if W))-/0*0) -M Z )}Pkd0Z(dx), (2.20) 

JeJo fceN \ i=i / 

with the convention that the sum over j is zero when k = 0. The generator £2 is given by: 

£ 2 <j>f( Z ) = (Z,Lf)^((ZJ)) + (Z,a(x)f /2 (x))r((Z,f)). (2.21) 

For a test function / : (t,x) 1— y ft(x) in C^' 2 (M-|- x M,1R), the evolution of (Zf, t > 0) can then be 
described by the following SDE: 

{z u f t )=uxl)+ f I i {ueVs _ } \J2fs(F^\x^,e))-f s (x^))p(ds,du,dk,de) 

JO i»xNx[0,l] \~{ I 

+ f f (Lf s (x)+d s f s (x))Z s (dx)ds+ f V V2a{X u s )d x f s {X u s )dB u s . (2.22) 
Jo Jr Jo u£Vg 

where (B u ) u ^u a family of independent standard Brownian motions. In [7J, such splitted diffusions 
are considered to describe a multi-level population. The cells, which correspond to the individuals in 
the present setting, undergo binary divisions, and contain a continuum of parasites that evolves as a 
Feller diffusion with drift b(x) = (b — d)x and diffusion <J 2 (x) = 2a 2 x. At the branching time s for the 
individual u, each daughter inherits a random fraction of the value of the mother. The daughters ul 
and u2 start respectively at f} 2) (X£ _,9) = G~ 1 (8)Xf_ and F 2 (2) (X S M _ , 9) = (1 - G~ 1 (9))Xf_, where 
G _1 is the generalized inverse of G, the cumulative distribution function of the random fraction. □ 



3 The auxiliary Markov process and Many- To-One formulas 

In this section, we are interested in the distribution of the path of an individual picked at random 
in the population at time t. By choosing uniformly among the individuals present at time t, we give 
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a more important weight to branches where there have been more divisions and more children since 
the proportion of the corresponding offspring will be higher. Our pathwise approach generalizes [T7j 
(discrete time) and [7] (continuous time Yule process). As mentioned in the introduction, this size bias 
has already been observed by [52^ [T2] for the tree structure when considering marginal distributions 
and by |14[ 155] for local branching Markov process. 

In Section 13. 1\ we introduce an auxiliary Markov process which approximates the distribution of 
an individual picked at random among all the individuals living at time t. The relation between Xj 
and the auxiliary process also appears when summing the contributions of all individuals of T (Section 
I3.2|) and of all pairs of individuals (Section 

3.1 Auxiliary process and Many- To-One formula at fixed time 

We focus on the law of an individual picked at random and show that it is obtained from an auxiliary 
Markov process. This auxiliary Markov process Y = (Y, A) has two components. The component Y 
describes the motion on the space E. The second component, A, encodes a virtual genealogy and 

Y can then be seen as the motion along a random lineage of this genealogy. More precisely, A is a 
compound Poisson process with rate rm ; its jump times provide the branching events of the chosen 
lineage and its jump sizes are related to the offspring number H, whose distribution is the size biased 
distribution of p. As for the motion, Y behaves like X between two jumps of A. At these jump times, 

Y starts from a new position given by Fj (., 0) where J is uniform on {1, . . . , H} and G is uniform 
on [0,1]. 

For the definition of A, we shall consider the logarithm of the offspring number as this is the 
quantity that is involved in the Girsanov formulas. Notice that we cannot recover all the jump times 
from A unless there is no offspring number equal to 1, that is p\ = 0. This can however always be 
achieved by changing the value of the jump rate r and adding the jumps related to F^ to the process 
X. Henceforth, we assume without loss of generality that: 

Assumption 3.1. The offspring distribution satisfies p\ = 0. 

By convention for a function / defined on an interval /, we set fc = (f(t),t G C) f° r an y C C /. 

Definition 3.2. Let Xj be as in Definition \2. 6\ with starting distribution [i G V(E). The corresponding 
auxiliary process Y = (Y,A), with Y = (Yt,t > 0) and A = (At,t > 0), is an E x M-valued cddldg 
Markov process. The process (Y,A) and I = (Ik,k G N*), a sequence of random variables, are defined 
as follows: 

(i) A is a compound Poisson process: At = J2k=i^°&(Hk), where S = (St,t > 0) is a Poisson 
process with intensity rm, and {H^,k G N*) are independent random variables independent of S 
and with common distribution the size biased distribution of p, (hph/m, h G N*). 

(ii) Conditionally on A, (Ik,k G N*) are independent random variables and Ik is uniform on 
{!,..., H k }. 

(in) Conditionally on (A, I), t\ = inf{t > 0; St ^ Sq} is known and the process YtQ^ is distributed 

aS X [0,T!)- 

(iv) Conditionally on (A,X, Yjo jTl )), Y Tl is distributed as F- (Y Tl -,0), where Q is an independent 
uniform random variable on [0, 1]. 
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(v) The distribution of (Y T1+ t,t > 0) conditionally on (A,I, Yro n ]) is equal to the distribution ofY 
conditionally on (A n+ t — A T1 ,t > 0) and (ii+fc, k G N*), and started at Y T1 . 

We write when we take the expectation with respect to (Y, A, I) and the starting measure is /U 
for the Y component. We also use the same convention as those described just after Definition 12.61 

The formula (|3.1|) in the next Proposition is similar to the so-called Many-to-One Theorem of 
Hardy and Harris [33] (Section 8.2) that enables expectation of sums over particles in the branching 
process to be calculated in terms of an expectation of an auxiliary process. Notice that in our setting 
an individual may have no offspring with positive probability (if po > 0) which is not the case in |33j. 

Proposition 3.3 (Many- To-One formula at fixed time). Fort > and for any nonnegative measurable 
function f G B(H([0,t],E xR),M + ) and t > 0, we have: 

jgF^j E^[f{Y m )\. (3.1) 

Remark 3.4. • Asymptotically, Nt and E[JVj] are of same order on {W > 0}, see (|2.1U|) . Thus, the 
left hand side of (|3.ip can be seen as an approximation, for large t, of the law of an individual 
picked at random in Vt- 

• For m > 1, a typical individual living at time t has prolific ancestors with shorter lives. For 
m < 1, a typical individual living at time t has still prolific ancestors but with longer lives. 

• If births are local (i.e. for all j < k, Fj k \x, 9) = x), then Y is distributed as X. 

Proof of Proposition [X3l Let A be a compound Poisson process as in Definition 13.21 (i). Let us show 
the following Girsanov formula, for any nonnegative measurable function g: 

E[g(A m )] =E[g(A! [QA )e^ m - 1 ^] , (3.2) 

where the process A' is a compound process with rate r for the underlying Poisson process and 
increments distributed as log(z^) with v distributed as -p. Indeed, g(A[ 0it ]) is a function of t, of the 
times T q = inf{t > 0; St = q} — inf{t > 0; St = q — 1} and of jump sizes log(H q ) of A: 

ff(A[ 0) t]) =^G g (i,ri,. . . ,r q ,Hi, . . . ,H q )l^q_ Ti <t<^?j' 1 T < }' 
for some functions (G q , geN). We deduce that: 
Eb(A [0 , t ])] 

+oo . q , 

/ E (rmye- rmt G q (t,t 1 ,...,t q ,h 1 ,...,h q )H^l {yq <t}< ft x . . . dt q 



q=0 + hi,...,h q i=l 



1 I. *± 

E / E r'e-^,^!,...,*^!,..^^-*- 1 )^'*)!]^^ k<t} d tl ...dt q 

q=0 jR9 + h lt ...,h q 8=1 ' 1 

E^Ajo^e-^- 1 ^ . 
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Recall that (St,t > 0) (resp. > 0)) is the underlying Poisson process of A (resp. A u ). Notice 

that if \u\ = q, then {S" = q} = {a(u) < t < /3(u)}. We thus deduce from (pT2|) that for q £ N, u G U 
such that \u\ = q, 



E[g(A m )l {St=q} ]=E g(Af 0jt] )e 



-r(m—l)t+A 



* 1 



{a(u)<t<0(u)} 



(3.3) 



Let q G N*. By construction, conditionally on {Ar 0) d = Ar ^i}, {5t = q}, . . . , I q ) = u}, Yi 0)t ] is 
distributed as -XTj^i conditionally on {Aj^ = A[ 0j t]}. This holds also for q = with the convention 
that (Ii, . . . , I q ) = 0. Therefore, we have for any nonnegative measurable functions g and /, 

E M [g(A [0jt] )/(Y [0j t])] = E 1 {|«|=g} E Mb(A[o,t])/(%t])l{(/ 1 ,...,/ 9 )=«}l{s t = g }] 



Using the points (i) and (ii) of Definition 13, 2\ we see that 

g 

P(h = u u ...,I q = u q \ A t , {S t = q}) = J] 1/H k = e' At if Ul < H u . . . u q < H q 

k=l 

= otherwise. 

Hence: 

E M [0(A [O)t] )/(Y [O)t] )] 

= E E l{|«|= 9 } E Mb(A [0 ,i])^[/(^ ] )|Af 



0,t]J|A« jt] =A [o>t] e ' 1 {«i<Hi,...,« 3 <H g } 1 {S t =g}] 



-r(m— l)t 



■{Mi<i/0,...,u q <i/ ( „ 1 .. )}l{5 t "=g}] 



thanks to (|3.3|) , Remark that since tt = «i . . . u g , {ui < u$, . . . , u q < V(u 1 ...u q -i)} = {i£ T}, and on 
this event, we have {5" = g} = {u E Vt} as noticed before. As a consequence: 

E M b(A [0it] )/(Y[ 0lt] )] = E E M^,t])f(Xio,t])e- rim - l)t 
Finally, we use a monotone class argument to conclude. ■ 



3.2 Many-to-Ones formulas over the whole tree 

In this section, we generalize identity (|3.1|) on the link between the tree indexed process Xj and the 
auxiliary Markov process Y by considering sums over the whole tree. 



Let us consider the space T> of nonnegative measurable functions / G 



.,Ex 



such that f(t,y) = f(t,z) as soon as y\o,t) = z [o.t)- By convention, if y is defined at least on [0, t), we 
will write f(t, 2/[o,t)) for f(t,z) where z is any function such that Z\q^ = U[o,t)- 

Proposition 3.5 (Many- To-One formula over the whole tree). For all nonnegative measurable func- 
tion f of V, we have: 



E,, 



))' 



ds e r( - m ~ l > E„ 



o 
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f(s,Y [0t 



(3.4) 



Before coming to the proof of Prop. 13.51 we introduce a notation that will be very useful in the 
sequel. By convention for two functions /, g defined respectively on two intervals If,I g , for [a, b) C If 
and [c, d) C I g , we define the concatenation [f[ a ,b) j 9[c,d)] =hj where J = [a, b + (d — c)), 



/(*) 



if t € [a, b) 



g{c+{t-b)) if t e [b,d- c + b). 



Proof of Prop. \3.b\ We first notice that if r is an exponential random variable with mean 1/r (r > 0), 
then we have, for any nonnegative measurable function g, 



E 



r / g(t)dt 
Jo 



E[g(r) 



(3.5) 



Besides, we have 



1 {u&T}f(l3{u),X^ u)) ) =E M l {ueT} f(/3{u), [^[o, a («));^[0,/3(«)-a(«))]) 



where conditionally on -Xp a ( u ))> /3("u), G 7"}, X = (X, c) with X of distribution Px u ( j and c 
the constant process equal to A£Vn. Notice that we have chosen X independent of (5{u). Thus, 
conditioning with respect to K a w);-^[o,+oo)]) { n £ 7~} and using (|3.5|) . we get 



E, 



We deduce: 



l {ueT}fW( u )> X [o,i3(u))) 



rE 



E, 



L {«eT} 



rE 



j3(u)—a(u) 



ds f{a{u) + s, X, 



[0,o(ii)+j)J 



l{«eT} / , ds f(s,X$ >a) ) 



a(u) 



Using Proposition 13.3} we get 



E„ 



52mu),x$ >m 



'rOO 

r I ds E 



hoo 

/• / ds E, 
lo 



l {ueV s }f(s,Xf^ s) ) 



u€V 3 



f 

JO 



-oo 



ds e r{m - 1)s E, 



The equality (|3.4|) means that adding the contributions over all the individuals in the Galton- 
Watson tree corresponds (at least for the first moment) to integrating the contribution of the auxiliary 
process over time with an exponential weight e r ( m-1 )* which is the average number of living individuals 
at time t. Notice the weight is increasing if the Galton- Watson tree is supercritical and decreasing if 
it is subcritical. 



Remark 3.6. We shall give two alternative formulas for 
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We deduce from (|3.4p that, for all nonnegative measurable function /, 



E„ 



£/(/3(u),Xf 0j/J(u)) ) 
.ueT 



E, 



f(r,Y [0:T) )e r 



(3.6) 



where r is an independent exponential random variable of mean 1/r. Thus, the right hand 
side of equation (|3.4|) can be read as the expectation of a functional of the process Y up to an 
independent exponential time r of mean 1/r, with a weight e rmT . 

Let T q = inf{t > 0; St = q} the time of the q-th jump for the compound Poisson process A. 
Using (|3.5|) . it is easy to check that, for any nonnegative measurable function g, 

1 ~ f+oo _ 

-£ E Mfo(%r,),T,)] =r / E4 5 (F [0)S) , S )] ds. 

771 g>l J ° 

Therefore, we deduce from (|3.4|) that, for all nonnegative measurable function /, 



E„ 



-E E m \fir q ,Y [0jTq) )e 



in 



9>1 



r(m-l)r q 



(3.7) 



This formula emphasizes that the jumps of the auxiliary process correspond to death times in 
the tree. 



3.3 Identities for forks 

In order to compute second moments, we shall need the distribution of two individuals picked at 
random in the whole population and which are not in the same lineage. As in the Many- To-One 
formula, it will involve the auxiliary process. 
First, we define the following sets of forks: 

TU = {{u,v) e U 2 : \u /\v\ < min(|u|,H)} and TT = TU n T 2 . (3.8) 
Let J2 be the operator defined for all nonnegative measurable function / from (E x M) 2 to M by: 

J 2 f(x,\)=[ 1 E Pkf(Fj l k Hx,0),X + log(k),Fj; k \x,0),X + log(k))de. (3.9) 

(a,fe)€(N*) 2 fc>max(a,fe) 

Informally, the functional J2 describes the starting positions of two siblings. Notice that we have 

J 2 f(x, A) = m (\[{H - l)/(if°(x, 6), A + log(iT), F { K H) (x, 8), A + \og(H))}d9, (3.10) 
J 

where H has the size-biased offspring distribution, and conditionally on H, (I, K) is distributed as a 
drawing of a couple without replacement among the integers {1, . . . , H}. 

For measurable real functions / and g on E x R, we denote by / <S> g the real measurable function 
on (E x R) 2 defined by: (/ (g> g)(x, y) = f(x)g(y) for x, y G E x M. 



15 



Proposition 3.7 (Many- To-One formula for forks over the whole tree). For all nonnegative measurable 
functions tp, ip G V, we have: 



E„ 



J 2 (e'. fo(t + t', [y m ;Y{ oy) ]) e— ' ] |4=t ~ =? 



(y T _) , (3.11) 



where, under E^ ; r is exponential with mean 1/r independent of Y , and, under ~E' X x , (Y',t') is 
distributed as ({Y, A + A), r) under K x . 

Proof. Notice that {(it, v) G Tli\ is equal to {3(u;, u, v) G W 3 ,3(a, 6) G (N*) 2 ,a ^ b,u = wau,v = 
wbv}. Let A be the l.h.s. of (13TTL) . We have: 

^=E E E + (/?(^) - /3HM^^ 

x ^(/?h + (p( wb v) - p( w )), [^(.));^ ) , /3( ^ )) ])l { ^ e r } 

Using the strong Markov property at time /3(w), the conditional independence between descendants 
and Proposition 13.51 we get: 



A=E E % 



E~ 



<p(i + r', [x [0it) ; ?jo,r')])' 



\t=P(w),x=X< 



x E'~ 



^(t + r', [5[ ,t); ?[o,t')]) 



l{«jaeT} 



|t=/8(ui),«=X' 



where under E^, A , (Y',t') is distributed as ((Y,A + A),r) under E x . As {wa,wb G T} = {w G 
T} fl {max{a, b} < u w } we have: 



A = E E M [l {w6T} J 2 (E'.[ ¥3 (t + r', [i^s^e"^] 



t=@(w),x=X v 



E'.ty(t + r / , [fc [0 , t) ;yj{ )iT , ) ]e r 



(3.13) 



with J 2 defined by (|3.9|) . The function under the expectation in (|3.13|) depends on [3(w) and X^^^. 
Equality (|3.6|) then gives the result. ■ 

We shall give a version of Proposition 13.71 when the functions of the path depend only on the 
terminal value of the path. We shall define J 2 a simpler version of J 2 (see Definition (|3.1U|) ) acting 
only on the spatial motion: for all nonnegative measurable functions / from E 2 to R, 



J 2 f(x)=m [\[(H-l)f(Fi H \x,e),F^\x,e))]d9, 
Jo 



(3.14) 



where (H,I,K) are as in (|3.10|) . 

The following Corollary is a direct consequence of Proposition 13.71 and the fact that Y is cadlag. 
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Corollary 3.8 (Many- To-One formula for forks over the whole tree). Let (Qt,t > 0) be the transition 
semi-group ofY. For all nonnegative measurable functions f,g € P, we have: 



E, 



[ ]T f(P(u),X^ u) _)g((3(v),XV (v) 



[ e r(m-l)(s + t + t>) ds dt dt , jlQs (J 2 (QJ g, Q vgv +s)) , (3.15) 

J[0,oo)3 



where ft(x) = f(t, x) and gt{x) = g(t, x) for t > and x £ E. 

We can also derive a Many- To-One formula for forks at fixed time. 

Proposition 3.9 (Many- To-One formula for forks at fixed time). Let t 6 M +; and ip,ip two nonneg- 
ative measurable functions on B([0, t],E). We have: 



M E ^[V])^(*[o,*]) 



r e 



2r(m-l)i 



-r(m-l)a da 



E, 



J 2 (E'. ^([y [0>a) ; ^, t _ o] ])] ,~ =? ® E'. ^([y [0>a ); ^, t _ a] ]] ,~ =? ) (^)] , (3.16) 



where, under K' x x , Y' is distributed as (Y, A + A) under E x . 

The l.h.s. of (|3.16|) approximates the distribution of a pair of individuals uniformly chosen from 
the population at time t. Indeed, we have in the r.h.s. of (|3.16|) an exponential weight eM™- 1 )* 
and thanks to Lemma 12.44 we know that K[Nt(Nt — 1)] ~ £7e 2r( - m ~ 1 ) i . The distribution of the paths 
associated with a random pair is described by the law of forks constituted of independent portions 
of the auxiliary process Y and splitted at a time a £ [0,t]. Notice that (|3.16p indicates that the fork 
splits at an exponential random time with mean l/r(m — 1), conditioned to be less than t. 



Proof of Proposition \3.9\ The proof is similar to the proof of Proposition 13.71 except that we use 
Proposition 13.31 instead of Proposition 13.51 to obtain an analogue of (|3.12|) . ■ 



4 Law of large numbers 

In this section, we are interested in averages over the population living at time t for large t. When 
the Galton- Watson tree is not supercritical we have almost sure extinction, and thus we assume here 
that m > 1. 



4.1 Results and comments 

Notice that Nt = implies Zt = and by convention we set Zt/Nt = in this case. For t S M + and 
/ a real function defined on E, we derive laws of large numbers for 

f(X?) , (Z t ,f) _Z»ev t f(X?) 
N N t E[N t ] E[N t ] ' l4 ' ij 

provided the auxiliary process introduced in the previous section satisfies some ergodic conditions. 
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Let (Qt,t > 0) be the semigroup of the auxiliary process Y from Definition 13.21 

Mf( Y t)\ = t*Qtf (4-2) 

for all /j, G V(E) and / nonnegative. Recall the operator J2 defined in (|3.14p . 

We shall consider the following ergodicity and integrability assumptions on /, a real measurable 
function defined on E, and \x G V(E). 

(HI) There exists a nonnegative finite measurable function g such that Qt\f\(%) < g(%) for all t > 
and x £ E. 

(H2) There exists 7r G V(E), such that (71", |/|) < +00 and for all x 6 E, Hm.t-t+00 Qtf(x) = (ir, f). 

(H3) There exists a < r(m — 1) and c\ > such that \1Qtf 2 < ci e a * for every i > 0. 

(H4) There exists a < r(m — 1) and C2 > such that [xQtJ2 S y9 ® 9) < C2 e Q * for every i > 0, with (7 
defined in (HI). 

Notice that in (H3-4), the constants a,c± and c% may depend on / and \x. 

Remark 4.1. When the auxiliary process Y is ergodic (i.e. Y converges in distribution to it £ V(E)), 
the class of continuous bounded functions satisfies (Hl-4) with g constant and a = 0. In some 
applications, one may have to consider polynomial growing functions. This is why we shall consider 
hypotheses (Hl-4) instead of the ergodic property in Theorem 14.21 or in Proposition 14.31 

The next Theorem states the law of large numbers: the asymptotic empirical measure is distributed 
as the stationary distribution ir of Y. 

Theorem 4.2. For any [i G V(E) and f a real measurable function defined on E satisfying (Hl-4), 
we have 

(ttJ)W mL 2 (P M ), (4.3) 
(n,f)l{W?0} in F ^-probability, (4.4) 

with W defined by h2.1b]) and it defined in (H2). 

For the proof which is postponed to Section 14.21 we use ideas developed in [T7j in a discrete time 
setting. We give an intuition of the result. According to Proposition 13.31 an individual chosen at 
random at time t is heuristically distributed as Yj, that is as tt for large t thanks to the ergodic 
property of Y (see (H2)). Moreover two individuals chosen at random among the living individuals 
at time t have a MRCA who died early, which implies that they behave almost independently. Since 
Lemma 12.41 implies that the number of individuals alive at time t grows to infinity on {W 7^ 0}, this 
yields the law of large numbers stated in Theorem 14.21 

Notice that Theorem 11.11 is a direct consequence of Theorem 14.21 and Remark 14.11 

We also present a law of large numbers when summing over the set of all individuals who died 
before time t. Recall that D t = J2u&T !{/?(«)<*} denotes its cardinal. 

Recall 5 in Definition [321 Notice that E[Sy = rmt. We shall consider a slightly stronger hypothesis 
than (H3): 

(H5) There exists a < r(m — 1) and C3 > such that K^[f 2 (Yt)St] < c^e at for every t > 0. 
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lim <#4 

E[N t ] 
lim <^> 

t-s-+oo N t 



Proposition 4.3. For any [i G V(E) and f a nonnegative measurable function defined on E satisfying 
(Hl-5), we have 

f( X B(u)-) 

km = (<K,f)W mL 2 P M , 4.5 

Euer/C^u)-) 1 ^)^} , A1 . ttj (AR , 

hm l { AT t>0 } p7 = (t,/)%/0} in ¥ ^-probability, (4.6) 

i— >+oo .L/£ 

defined by \2.1(J\) and tt defined in (H2). 



We can then extend these results to path dependent functions. In particular, the next theorem 
describes the asymptotic distribution of the motion and lineage of an individual taken at random in 
the tree. In order to avoid a set of complicated hypothesis we shall assume that Y is ergodic with 
limit distribution tt and consider bounded functions. 

Theorem 4.4. We assume that there exists tt G V(E) such that for all x G E, and all real-valued 
bounded measurable function f defined on E, \\m.t^ f00 Qtf(x) = (TT,f). 

Let T > 0. For any real bounded measurable function ip on D([0, T],E x R+), we have 

1 



felfi^£^-T, t ],AjU,4^ W inL 2 (F,), 

1 i m W £ ^( X [t-T,t] > A ft-T,t] - A t-T) = E tt [</?(?[o,T])l l{wvo} in F ^-probability, 



1 u&Vt 

with W defined by QEQ) . 

Let J\ be the following operator associated with the possible jumps of Y: for all nonnegative 
measurable function / from E to R, 

J 1 f(x)=m {\[f(Ff\ X ,e))}d6, (4.7) 
Jo 

where H has the size-biased offspring distribution, and conditionally on H, I is uniform on {1, ... , H}. 

Proposition 4.5. We assume that there exists tt G V(E) such that for all x G E, and all real-valued 
bounded measurable function f defined on E, lim^oo Qtf(x) = {tt, /). 

Let ip be a real bounded measurable function defined on E-valued paths. We set, for x G E, 
f{x) = Ea;[</?(Y[o,Ti))]> with 7~i from Dehnition \S.2\ We have 



[a(«),/3(M))) 1 {/3H<0 



hm -!vr-r>\ = { 7r ,Jif)W inL ir^), (4i 



lim l{Ar t>0 } £ = (vr, Jif)l{wM\ in probability, 

m'tfi W defined by l2AU\) . 



Remark 4.6. The hypothesis on Y in Theorem 14.41 and Proposition 14.51 is slightly stronger than the 
ergodic condition (i.e. Y converges in distribution to tt), but it is fulfilled if Y converges to tt for the 
distance in total variation (i.e. for all x G E, linn._j.oo sup^g^- |P x (i* G A) — tt(A)\ = 0). This property 
is very common for ergodic processes. 
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4.2 Proofs 

Proof of Theorem \4.ty We assume (Hl-4). We shall first prove (|4.3p for / such that (it, f) = 0. We 
have 

■(ZtJ) 2 - 



E„ 



At + Bt 



where 



A t = E[N t ]- 2 E,[Y, f\X?)} and B t = E[N t ]- 2 E,[ ]T f(X?)f(X? 

uev t {u,v)ev t 2 



Notice that 



A 



t = e 



-r(m-l)t ]g 



-r(m—l)t 



vQtf 



t— too 



+ 0, 



(4- 



thanks to ()2.8|) and (|3.ip for the first equality and (H3) for the convergence. We focus now on B%. 
Notice that Proposition 13.91 and then (HI) and (H4) imply that 



E[N t ]- 2 E,[ ]T =r[ i*QsJ2(Qt-.\f\®Q±-s\f\)e- r{mr - 1)a ds 

is finite. We thus deduce that 



Bt = rj t fiQsJi (Qt-sf ® Qt-sf) e~^ m ~ l > ds. 



(4.10) 



Now, since (tt, /) = 0, we deduce from (H2) that for s fixed, and y,z £ E, lim i ^ 00 (Q 4 _ s / (g) 
Qt- S f)(y, z) = 0. Thanks to (HI), there exists g such that l^ s<t j\(Qt- s f <8> Qt_ s /)| < (g ® g) and 
(H4) implies that / °° ds e~ r ^ m ~ 1 ^ s fiQ s J2(g <8> 5) is finite. Lebesgue Theorem entails that 

lim £ t = lim r [ fiQ s J 2 (Qt-sf ® Qt- S /) e^™" 1 ^ = 0. 

This ends the proof of (|4.3|) when (tt, f) = 0. 
In the general case, we have 



{ZtJ) (*J)w= {ZtJ -^ f)) + (.,f) 



E\N t 



E[N^ 



E[N t _ 



W 



(4.11) 



Notice that if / and fi satisfy (Hl-4) then so do / — (tt, f) and fi. The first term of the sum in the 
r.h.s. of (|4.1ip converges to in L 2 thanks to the first part of the proof. The second term converges 
to in L 2 thanks to Lemma 12.41 Hence we get (|4.3p if / and fi satisfy (Hl-4). 

We deduce (031) from JO} and (I2TTU|) . ■ 

Proof of Proposition We assume (Hl-5). We shall first prove (|4.5p for / such that (71", f) = 0. 
We have: 



E\DA- 2 E, 



u) — l - 1 -{/3(«)<t} 



A t + B t + C t 
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where 



B t = E[Ar 2 E, 



C t = 2E[A]~ 2 E 



u) — > {P(u 



)<*} 



.ueT 



v)—l {P(u)<t,P{v)<t) 



v)—) *-{/3(v)<t} 



The terms At and Bt will be handled similarly as in the proof of Proposition 14.21 Notice that 

At = rE[D t r 2 ^dse^>E,[f(Y s ^] = £ ds^» pQ.f ^ 0, (4.12) 



thanks to (|3,4|) for the first equality, (|2.12|) for the second and (H3) for the convergence. 
Notice that Corollary 13.81 and then (HI) and (H4) imply that 



E 



/' X \f( X ^u)-)\\f{ X l(v)-)\ 1 {l3(u)<tfi(v)<t} 



'[0,+oo) 3 

is finite. We thus deduce that 
r 3 (m — l) 2 



3 M<3 S ^(Q S '|/I ® Qs"\f\) e 



r(m-l)(s+s'+s") 



l{ s+s > <t , s +s"<t} dsds'ds' 1 



B 



1 " {e{m-i)t _ 1)2 j %+oo) , 

r 3(jn — I) 2 e 2r(m— l)t 



MQs-4(Q s '/ ® Q s "/J e 



r(m-l)(s+s'+s") 



l{ s+s /< tjS+s //< t } dsds'ds' 1 



(r ,,,„_ J , _ , )2 y +oo)3 »QsJ2(Qt-t>f ® <?*-*»/) e^— 1 )^-*'-*") l W<M<t » <i} tfadt'dt". 

Now, since (vr, /) = 0, we deduce from (H2) that for t',t" fixed and y, 2: G £7, limi_ ) . 00 ((5t_i// (g) 
Qt-t"f)(y, z ) = 0. Thanks to (HI), there exists g such that | (Qt—f/ <£> Qt-t"f)\ — (<7® <?)• Then (H4) 
implies that Jj +00 )3 /xQ s J2(g <8> 5) e^" 1-1 ^ - * ~* ) l{ s <t', s <t"} dsdt'dt" is finite. Lebesgue Theorem 
entails that 

lim 5t = 0. (4.13) 

t— >oo 

Let us now consider Ct ■ We have Ct < C{ + C" where 



Ct = E[A] _2 E At 



X / 2 ( X ^)~) 1 {/3M<0 



and C" = E[A] E 



X / 2 W u)—) *-{P(v)<t} 



We deduce from (13.41) that 



c; = e[a]~ 2 e, 



ElH/ 2 (^)-)i w , 



.»eT 



E[A]~ r / e' 
Jo 



r(m— l)s 



E„ 



)<*} 



S s -f 2 (Ys-] 



( e f(m-i)*_i)2 y 



5 s / 2 (y a ) 
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We deduce from (H5) that 



lim C[ = 0. 



t— >oo 



(4.14) 



Using the conditional expectation w.r.t. X u , (|2,12p and (|3.4|) . we get 



C'l = E[A]" 2 E M 



E / 2 ( X |(«)~) 1 {/3(«)<O mE E 1 {^(«)<t'}]|*'=*-^(u) 



m — 1 
rn 



E ^ 2 (^(«)-) 1 {/3(«)<o( e 



r(m-l)(t-/3(u)) 



-I) 



E[D t ]~ 2 / ds e r(m - 1)s E„ / 2 (y s _)(e r(m - 1)(i - s) -1) 
m-l J L 



m(m - 1) e^™" 1 )* f* „ „ 9 



( e r(m-l)t _!)2 7 Q 

We deduce from (H3) (or (H5)) that 



lim C'l = 0. 



t— >oo 



(4.15) 



The proof of (|Q]> . when (vr, /) = 0, is then a consequence of (|4TT2l . (l4TTTfl) . (14TT4T) and (JUSJ. 
In the general case, we have 



eia]- 1 E /(^( U )-) i {^m<*} - (*,f)w 



nDt}- 1 E (/(*««)-) - («•. /)) iw-X*} + /) (j^j - W) . (4.16) 



Notice that if / and [i satisfy (Hl-5) then so do / — {tt, f) and [i. The first term of the sum in the 
r.h.s. of (|4.16|) converges to in L? thanks to the first part of the proof. The second term converges 
to in L 2 thanks to Lemma 12.51 Hence we get (|4.5p if / and \x satisfy (Hl-5). The convergence in 
probability is thus obtained thanks to (|4.5|) and (|2.10|) . ■ 



Proof of Theorem \4-4\ The proof is similar to the proof of Theorem 14.21 Some arguments are shorter 
as we assume that <p is bounded. 

We shall first consider the case E T [(^(Y[o,T])] = 0- We assume that t > T. We have 



u&Vt 



A t + B' t + Bl 



where 



A t = E[N t ]- 2 E,[ E V 2 (X?t-T,t], ^t-T,t] ~ K-t) 



ueVt 



(u,v)€V t 2 



B'l = E[Ar 4 r 2 E M [ E <P(*lt-T,t], A ft-T,t] " K-TMXf t -T,t]^[t-T,t] ~ K-t)1 



{0(uAv)>t~T} 



{0(uAv)<t-T} 



(u,v)€V t 2 
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We assume that ip is bounded by a constant, say c. We have At < c 2 E[Nt] 1 so that lim^oo At = 0. 
We have, using Proposition 13.91 



1 5; | <c 2 E[^]-%[ J2 

L {/3(uAv)>t-T} 

(u,v)eV t 2 

so that limt^oo B[ = 0. 

We set /(a?) = E x .[(/9(Y[ 0i t])]. Using Proposition 13.91 once more, we get 

B't' = r f e-K-D« E, [j 2 ( E '. ^(^_ a _ T , t _ a] , A| t _ a _ T)t _ a] - A' t _ a _ T )] 



IE'. [v(^«-T,t-o]. A| t _ a _ T)t -a] " K-a-T)]){Ya 



t-T 

r I e 
/o 



-r(m— l)a 



HQaJ2(Qt-a-Tf <8> Qt-a-Tf)da 



l{a<t-T}da 



By hypothesis on y, we have that, for fixed a, lim^oo Q t - a -Tf = (vr, /) = 0. Using Lebesgue 
Theorem, we get lim^oo B" = 0. This gives the result for the L 2 (P^) convergence when (vr,/) = 0. 
We conclude in the general case and for the convergence in probability as in the proof of Theorem 



Proof of Proposition \4-5[ The proof is similar to the proof of Proposition 14.31 Some arguments are 
shorter as we assume that ip is bounded. 

We shall first prove (|4.8|) for <p such that (tt, J\f) = 0. We have: 



E[Ar'E 



where 



A t = E[A]" 2 E M 



B t = E[A1 E 



t P( X [a(u)Mu))) 1 {P(u)<t} ) =A t + B t + C t , 

,«eT / 



V(^(«),/9(«))) 2l {i9(«)<t} 



^2 V 3 (^(«),/3(u)))v ::> (^(t ) ),/3(,;))) 1 {^(«)<t,/3(D)<t} 

(«,^)eJ r T 



C t = 2E[A]-% 



H ( P( X [a(u),(5(u))) l f( X [a(v),l3(v))) 1 Wv)<t} 



We assume that 99 is bounded by a constant, say c. We have A t < c* fE[D t ] so that lim^oo A t = 0. 
Thanks to Corollary 13.81 we have 



\C t \ < 2c 2 E[D t ]- 2 E^[J2 \v\l{j}(v)<t}] 

= 2c 2 E[D t }- 2 f ds e^-^EJSs] 
Jo 

= 2c 2 E[D t }- 2 f ds srme r(m - 1)s . 
Jo 
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This implies that lim^oo Ct = 0. 

We set ht(x) = K x [(p(X< 0iT \)lr T<t \], where r is an exponential random variable with mean 1, 
independent of X. 

Using the conditional expectation w.r.t. X u , where v! is the ancestor of u, and X v , where v' is 
the ancestor of v, we have, according to u' = v' or u' ^ v', 



B t = B' t + B", 



where 



B' t = E[D t ]- 2 K, 



^2 Jz( h t-I3(u>) ® ^i-/3(«'))(^l(«')-) 1 {/ 3 ("')<*} 



u'eT 



B'l = E[A]~% 



Y2 ^i(^-/3(«o)(^(«o-) Ji (^-^(«o)^(v')-) 1 {) 9 ( tt ')< / " ' 'i 

Using the definition of J 2 , (pOlj) . we get < c 2 E[A] _1 (? 2 + 

m 2 — m) and thus lim£_>. o0 = 0. 

We deduce from Corollary 13, S\ that 
r 3 (m — l) 2 



B'l 



dsds ds ' 



( e r(m-l)t_ 1 )2 7 [0 +oo)3 

nQshi^Qs'Jiht-s-s' ® Q s "Jiht~s-s") e r ( m-1 ^ s+s +s ^ 1{ S +,' /..s+.«". /} 
r 3 (m-l) 2 e 2r ( m - 1 ) < 



dsdv'dv' 1 



( e r(m-l)t _1)2 7 [0i+oo)3 

nQsJ2(Qt-s-v'Jih v ' ® Qt-s-v" J\h v 'i \ e - r ( m ~ 1 )( s + v + v ) l{ v / <t _ S)U » <t _ s }.. 

By hypothesis on Y, we have that, for fixed s and v, lim^oo Qt-s-vhh v = (vr, J\h v ). Using Lebesgue 
Theorem, we get 



lim B'l = r 3 (m - l) 2 
t— >oo 



[0,+oo) 3 



dsdv'dv" {TT,J 1 h v/ ){TT,J 1 h v// )e- r{m - 1){s+V ' +V " ) 



Notice that h t (x) = ±E x [ip(Y [0jTl) ) e^ 1 ^ l {n<t} ] so that 

/■+oo 1 

r(m - 1) ^ h t (x) e-^' 1 ^ = - E x [<p(Y [0tTl) )]. 
Recall f[x) = K x [ip(Y [0jTl) )]. We get lim^oo B'( = {m \ )m i (tt, Ji/) 2 = 0. Therefore, we get that 



lim A t + Bt + C t = 0, 



t— >OD 



which gives the result for the L 2 (P^) convergence when (tt, J\f) = 0. We conclude in the general case 
and for the convergence in probability as in the proof of Proposition 14.31 ■ 



5 Examples 

We now investigate several examples. In Section 15.11 splitted diffusions are considered as scholar 
examples. In subsection 15. 2[ we give a biological application to "cellular aging" when cells divide in 
continuous time, which is one of the motivation of this work. In Section 15.31 we & ve a central limit 
theorem for nonlocal branching Levy processes. 
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5.1 Splitted real diffusions 

A first example consists in binary branching: the continuous tree T is a Yule tree. For the Markov 
process X, we consider a real diffusion with generator: 

Lf[x) = b(x)f'(x) + ^f'(x). (5.1) 

We assume that b and a are such that there exists a unique strong solution to the corresponding SDE, 
see for instance [2S] Theorem 3.2 p. 182. 

When a branching occurs, each daughter inherits a random fraction of the value of the mother: 

F«(x,0) = G- 1 (9)x, F®(x,9) = (1 - G~ 1 (6))x, 

where G is the cumulative distribution function of the random fraction in [0, 1] associated with the 
branching event. We assume the distribution of the random fraction is symmetric: G{x) = l—G(l—x). 
The infinitesimal generator of Y is characterized for / £ C^(M,1R) by: 

Af(x) =b{x)f'{x) + a{x)f"{x) + 2r £ (^{f (G' 1 (9)x) - f(x)) + i(/((l - G^ 1 (9))x) - /(x))) d9 

=b(x)f'(x) + a(x)f"(x) + 2r / (f(qx) - f(x))G(dq). (5.2) 

J o 

Particular choices for the functions b and a are the following ones: 

(i) If b(x) = and <r(x) = a, we obtain the splitted Brownian process. 

(ii) If b(x) = —/3(x — a) and o~(x) = a, we obtain the splitted Ornstein-Uhlenbeck process. 

(iii) If b(x) = 1 and cr(x) = 0, the deterministic process X can represent the linear growth of some 
biological content of the cell (nutriments, proteins, parasites...) which is shared randomly in the 
two daughter cells when the cell divides. More precisely here, each daughter inherits random 
fraction of this biological content. 

Let us note that if b(x) = fix and a(x) 2 = a 2 x, we obtain the splitted Feller branching diffusion. But 
in this case, almost surely, the auxiliary process either becomes extinct or goes to infinity as t — > oo. 
The assumption (H2) is not satisfied. This process is studied in [7] as a model for parasite infection. 

The following results give the asymptotic limit of the splitted diffusion under some condition which 
is satisfied by the examples (i— iii) . For this we use results due to Meyn and Tweedie [45} |4"6] . 

Proposition 5.1. Assume that Y is Feller and irreducible (see J^6] / p. 520) and that there exists 
K G M+, such that for every |x| > K, b{x)/x < r' with r' < r. Then, the auxiliary process Y with 
generator A is ergodic with stationary probability ir. Furthermore Y^ u£ y t 5x^(dx)/Nt converges weakly 
to tt as t — > oo and this convergence holds in probability. 

Proof. Once we check that Y is ergodic, then Corollary 11.21 and the fact that W defined by (|2.10p is 
a.s. positive readily imply the weak convergence of the Proposition. To prove the ergodicity of Y, 
we use Theorems 4.1 of [45J and 6.1 of |46j. Since Y is Feller and irreducible, the process Y admits 
a unique invariant probability measure n and is exponentially ergodic provided the condition (CD3) 
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in is satisfied, namely, if there exists a positive measurable function V : x i— > V(x) such that 
lira a ._ + ± oc V(x) = +00 and for which: 

3c> 0, d G R, Vx G R, AF(x) < -cV(s) + d. (5.3) 

For V(x) = |x| regularized on an e-neighborhood of (0 < e < 1), we have: 

V|x| > e, AV(x) =sign(x)6(x) + 2r\x\ / (q — l)G(dq) = sign(x)6(x) — r\x\, (5-4) 

Jo 

as the distribution of G is symmetric. By assumption, there exists 77 > and K > e, such that f|5.4|) 
implies: 

Vx G R, AV(i) < -r]V(x) + ( sup |6(a?)| + rK)ln xl < K} . (5.5) 

This implies (|5.3p and finishes the proof ; the geometric ergodicity expresses here as: 

3/3 > 0, B < +00, Vt G R+, Vx G R, sup |Qt5(a;) - (vr, £?> | < 5(1 + |x|) e _/3i . (5.6) 

9/|s(tOI<i+M 

■ 

Remark 5.2. The examples (i-iii) satisfy the assumptions of Proposition 15.11 If b and a are bounded 
Lipschitz functions, X is Feller (e.g. Theorem 6.3.4 p. 152 of [54] ) . and thus Y is also Feller. The 
Feller property also holds for Ornstein-Uhlenbeck processes. The irreducibility property is well known 
for diffusions as (i) and (ii) and trivial for (iii). 

Remark 5.3. If there exists K > in Prop. 15.11 such that for every \x\ > K, 2b(x)/x + Qa(x)/x 2 < r' 
with r' < rj \l - q 4 ) 2 G(dq), then we can use similar arguments as in the proof of Proposition 15.11 
We get that the auxiliary process Y is geometrically ergodic with 

3/3 > 0, B < +00, Vt G R+, Vx G M, sup \Q t g(x) - (vr, g) \ < B{1 + |x| 4 ) e - ^ (5.7) 

Wl9(«)l<i+M 4 

instead of (|5.6p . This result will be used for the proof of the central limit theorem. 
5.2 Cellular-aging process 

We now present a generalization to the continuous time of Guyon [32] and Delmas and Marsalle [T7j 
about cellular aging. When a rod shaped cell divides, it produces a new end per progeny cell. So each 
new cell has a pole (or end) which is new and an other one which was created one or more generations 
ago. This number of generations is the age of the cell. Since each cell has a new pole and an older one, 
at the next division one of the two daughters will inherit the new pole and the other one will inherit 
the older pole. Experiments indicate that the first one has a larger growth rate than the second one 
(see Stewart et al. [SB] for details), which indicates aging. 

To detect this aging effect, [321 E] used discrete time Markov models by looking at cells of a 
given generation. To take into account the asynchrony of cell divisions, it may be useful to consider 
a continuous time genealogical tree. 

We consider the following model. Cells are characterized by a type 77 G {0, 1} (type corresponds 
to a cell of age 1 and type 1 to cell of greater age) and a quantity ( (growth rate, quantity of damage 
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in the cell) that evolves according to a Markov process depending on the type of the cell. Cells may 
die, which leads us to the following model. At rate r, each cell is replaced by one cell of type (resp. 
1) with probability po >0 (resp. p\ > 0), by two cells of type and 1 with probability po,i > 0, or by 
no cell with probability 1 — po — p\ — po,i > 0. The way the quantity Q is given to a daughter depends 
on its type and on the fact that it has or has not a sister. 

This can be stated in the framework of Section 12.21 and Section [3l For the sake of simplicity, we 
shall assume that C, evolves as a real diffusion between two branching times. 

Let L° and L 1 be two diffusion generators: for / G C 2 (M. x {0, 1},M): 

L v f((,v) = b(C,v)d(f((,7l)+v(C,v)di ( f((,7l), V 6 {0,1}. (5.8) 

We assume there exists a unique strong solution to the corresponding two SDE, see for instance [36J 
Theorem 3.2 p. 182. We consider the underlying process X = ((Ct,Vt),t > 0) with generator 

Lf((, V ) = l {v=0} L°f((, 0) + 1 {??=1} L 1 /(C, !)• 

Notice the process (rjt,t > 0) is constant between two branching times. The offspring distribution is 

p(dk) = (1 -p - pi - Po,i)So(dk) + (po +pi)S 1 (dk) + po,i5 2 (dk). (5.9) 

The offspring position is given by: 

if } ((C, 77), 9) =(<7o(C, V ,9 X ), 0)l { e 2 < P0/iP0+Pl)} + (gi (C, V,0x), 1)1 { * 2>Po/(PO+Pl)} 
F8?l((C,l),0)=(9?\<:,V,O),i) for i€ {0,1}, (5.10) 

(2) (2) 

for some functions go,gi,gQ ,g\ and (6\ , 62) a function of 9 such that if 9 is uniform on [0, 1] , then 9\ 

(2) (2) 

and #2 are independent and uniform on [0, 1]. The division is asymmetric if ^ g\ . One important 
issue is, using the law of large number (Section H]) and fluctuation results, to test if the division is 
asymmetric, which means aging, or not. Let us mention that a natural question would be to give the 
test in a more general model in which the division rate depends on the state of the cell and of the 
quantity of interest £ (which is realistic if for example £ describe the quantity of damage of the cell). 

Let us consider a test function / : (t, (,rj) (->■ ft(C,v) in ^' 2 (M + x(Ix {0, 1}),R), and let (B u ) u< z U 
be a family of independent standard Brownian motions. The SDE describing the evolution of the 
population of cells then becomes with the notations of (|5.8|) . (|5.9|) and (|5.1U|) : 



(Z t ,f t ) = (Z ,f )+ f [ l {ueVs } [J2fs(F^\(Ct,rit),0)) - f s (C_)) P(ds,du,dk,d9) 

JO iWx{0,l,2}x[0,l] \~{ I 

+ f I (L r >f(C,v) + dsfs(C,v))Z s (d(,dr ] )ds (5.11) 

JO itx{0,l} 

+ f E \jMQ,Vt)dcf(C,V^)dBl (5.12) 



If Y is ergodic, then l/jv t >o} J2 U £V t $Xf(dx)/N t converges to a deterministic non degenerated mea- 
sure on M + x {0, 1}. Given a particular choice for the parameters 

90, 9i, 9o\ 9i \ L°, L\ Po, pi and 
Po,i of the model, one can use arguments similar to the ones used in Proposition 15.11 and Remark 15.21 
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to prove the ergodicity of Y. Let us give examples. 



The first illustration that can be viewed as a direct generalization of the model of Delmas and 
Marsalle [T7]. The quantity £ £ R models the cell growth rate, which is assumed constant during the 

(2) (2) 

cell's life: b(Q,rj) = and cr(C,r]) = 0. For the functions go, g\, g$ , g\ which describe the daugthers' 
growth rates, as functions of their mothers' characteristics, we set: 

ffo(C,?7,#) = «oC + A) + £o, 9i{C,,ri,9) = a>i( + Pi + e 1 

g ( o ) (C,V,0) = a' C + (3' + e' , gf\^e) = «iC + # + 4, (5-13) 

where ao> c*i> o^, £ (—1, 1) and where /3o, /3 0) /^i £ R- The random variables Eq, £i> and e'i 
generated thanks to the uniform variable 6 and their distributions is as follows: e'q and e'-y are Gaussian 
centered r.v. with variances o\ > and cr 2 > respectively, while (eo,ei) is a vector of Gaussian 
centered r.v. with covariance 

a 2 >0, p£(-l,l). 

In [17], this model is used to test aging phenomena, for instance, which correspond to (ao,/3o) = 
Delmas and Marsalle in discrete time prove that the auxiliary process, which correspond 
to the Markov chain associated here to the continuous time pure jump process Y, is ergodic. As a 
consequence, Y is recurrent, admits an invariant probability distribution (since the jump rate r(po + 
pi + 2po,i) is a constant) and is hence ergodic (see e.g. Norris [48]). 

5.3 Branching Levy process 

We consider particles moving independently on R following a Levy process X and reproducing with 
constant rate r. Each child jumps from the location of the mother when the branching occurs. We 
are interested in the rescaled population location at large time. 
The generator of the underlying process X is given by: 

Lf{x) = bf'(x) + ^-f"(x) + / {f{x + y)- f{x) - yf'(x)l M<1} )h(dy) (5.14) 
2 Jr\{o} 

with b £ R, <7 £ R + and h a measure on R \ {0} such that / K y{ } V 2 h(dy) < +oo. The particles 
reproduce at rate r in a random number of offspring distributed as p = (pt,k £ N), such that 
J2k>i kPk > 1 (supercritical case). The offspring position is defined as follows: 

Fj k \x,e)=x + A*{9), je{l,...,k}, (5.15) 

where we recall that x is the location just before branching time and k is the number of offspring. We 
assume the following second moment condition: X^fceN^fc Sj=i E[A* (0) 2 ] < oo, where is uniform 
on [0,1]. 

Proposition 5.4. We have the following weak convergence in Mp(M): 

lim -!- 8 x u_p t {dx) = irz(dx)l( W>0 }, in probability (5.16) 
t— ¥+oo rs-i *■ — ' 1 _ — L ' 

* ueVt ^ 
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where 7T£ is the centered Gaussian probability measure with variance E and 

+00 k 



(3 = b+ yl {lym h(dy)+rJ2Pk^2n^m, (5-17) 

JR\{0} k=1 j=1 

„ +00 k 

/ y 2 h(dy)+rY,PkJ2 E ^j(®n (5-18) 



/K\{0} fe=1 j=1 

Proof. The auxiliary process Y is a Levy process with generator: 



a 2 



Af(x) = bf'(x) + —f"{x) + / {f{x + y)- f{x) - yf{x)l M<l} ) h{dy) 

fe=l m J o j=l fc 

In particular, we have for all iEK: 

„ +00 k 

E x [Y t ] = x + t(b+ yl wm h(dy)+rJ2Pk 1 £ E l A j( e ^)= x + l 3t ( 5 ' 19 ) 

JR\{o} k=l j=1 

„ +00 k 

E x [Y t 2 ]-E x [Y t ] 2 =t(<r 2 + / v i h(dy)+r y £p k y £M[A^e) 2 ]) = St. (5.20) 



fc=i j=i 

Then, we deduce from the central limit theorem for Levy processes or directly from Levy Khintchine 
formula, that ((Y t - /3t)/y/t,t > 0) converges in distribution to 7T£. This implies that for any fixed s, 
{(Y t - S — (3t)/y/i,t > 0) converges in distribution to 7T£. 
Let <p be a continuous bounded real function and define 

f t {x) := <p({x - Pt)/Vt\ fort>0, i£i 

Let (Qt, t > 0) be the transition semi-group of Y. We get that for any fixed s and 

Jim Qtsft(x) = {irs,(p}. (5.21) 

It is then very easy to adapt the proof of Theorem 14.21 with / replaced by ft — (vrs , <p) : (|4.9|) holds 
since is uniformly bounded; (|4.1U|) holds using similar arguments with (|5.21|) instead of (H2) and ft 
uniformly bounded instead of (HI) and (H4) arguments. Similar arguments as in the end of the proof 
of Theorem 14.21 imply that for any continuous bounded real function ip, the following convergence in 
probability holds: 

uGVt 

This gives ([536]) . ■ 
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6 Central limit theorem 



6.1 Fluctuation process 

In order to study the fluctuations associated to the LLNs, Theorem 14.21 we shall use the martingale 
associated to Z t , see (|2.17p . We focus on the simple case of splitted diffusions developed in Section 
15. 11 Our main result for this section is stated as Proposition 16.41 

In the sequel, C denotes a constant that may change from line to line. We work in the framework 
of Section 15.11 



We consider the following family, indexed by T > 0, of fluctuation processes. For / € 
and t > 

v^^-i^) (6,) 

where we recall that Nt = Card(Vj) = (Zt, 1} and Qt has been defined in (|4.2p . The family Qt is the 
transition semigroup of the auxiliary process Y, which is given by: 

Y t = X + fb{Y s )ds + [ a(Y s )dB s - f (1 - q)Y s _p(ds,dq) (6.2) 
JO Jo Jo 

where Xq is an initial condition with distribution fi, where (Bt)teR+ is a standard real Brownian 

motion and where p(ds,dq) is a Poisson point measure with intensity 2r ds (g) G{dq) with G such that 

Jj q f(q)G(dq) = Jj 1 j(/(g)/2 + /(l — q)/2)G(dq). As in Section |5TT| we will assume in the sequel that 

G is symmetric. In this case, G(dq) = G(dq). 

The idea in (|6.ip is to compare the independent trees that have grown from the particles of Zt between 
times T and t + T, with the positions of independent auxiliary processes at time t and started at the 
positions Zt- We recall that L is the generator defined in (|5.ip . and let J be the operator defined on 
the space of locally integrable functions by 

Jf{x)=--f(x)+r J q (f(qx) + f((l-q)x))G(dq) = - T f{x) + 2r J f{qx)G(dq). (6.3) 

This operator will naturally appear when computing the equation satisfied by rj T by applying (|2.22|) 
with f t (x) = e- rt / 2 f(x): 

Proposition 6.1. The fluctuation process &6.1\) satisfies the following evolution equation: 

(vLf) = f I (Lf(x) + Jf(x)) V ?(dx)ds + M?(f), (6.4) 
where Mf(f) is a square integrable martingale with quadratic variation: 



r I (f(qx) + /((l - q)x) - f{xjfG{dq) + 2a 2 {x)f'{xf 



^[N s+T \ 1 ' ■ 1 

The proof of this proposition is given in Section 16.31 In the following, we are interested in the 
behavior of the fluctuation process when T — > +oo. The processes r] T take their values in the space 
A4s(M) of signed measures. Since this space endowed with the topology of weak convergence is 
not metrizable, we follow the approach of Metivier [U] and Meleard [33j (see also |25t I55j ) and 
embed A^s(M) in weighted distribution spaces. This is described in the sequel. We then prove the 
convergence of the fluctuation processes to a distribution-valued diffusion driven by a Gaussian white 
noise (Proposition 16. 4p . 
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6.2 Convergence of the fluctuation process: the Central Limit Theorem 



Let us introduce the Sobolev spaces that we will use (see e.g. Adams [1]). We follow in this the steps 
of [H1I13]. To obtain estimates of our fluctuation processes, the following additional regularities for b 
and a are required, as well as assumptions on our auxiliary process. 

Assumption 6.2. We assume that: 

(i) b and a are in C 8 (R, R) with bounded derivatives. 

(ii) There exists K > such that for every \x\ > K, 2b(x)/x + 6o~(x)/x 2 < r' with r' < r J (l — 
q 4 ) 2 G(dq). 

(Hi) Y is ergodic with stationary measure tt such that (ir, \x\ 8 ) < +oo. 

(iv) for every initial condition \i such that (/i, |x| 8 ) < +00, sup 4gK+ E^JYj 8 ] < +00. 

Remark 6.3. (i) Notice that under Assumption 16.21 (i), there exist b and a > s.t. forall 1 6 I, we 
have \b(x)\ < 6(1 + \x\) and \a(x)\ < a(l + \x\). 

(ii) Conditions for the ergodicity of Y have been provided in Proposition 15.11 and Remarks 15.21 and 
5.31 Under Assumption 16.21 (ii) , Remark 15.31 applies and we have geometrical ergodicity with (|5.7|) . 

(iii) The moment hypothesis of Assumption 16.21 (iv) is fulfilled for the examples (i-iii) of Section 15.11 
provided the initial condition satisfies (/u, |x| 8 ) < +00. This can be seen by using Ito's formula (e.g. 
[55] . Th. 5.1 p. 67) and GronwalPs Lemma. Moreover, for every p £ {1, . . . ,8}, E^jjYtp'] < +00. 

(iv) Assumption 16.21 (iii) and (iv) imply: Vp £ {1, . . . , 7}, f R \x\ p ir(dx) < +00 and lim t _>. +00 E^[|Y| P ] = 
Ir \x\ p ir(dx). This is a consequence of the equi-integrability of (|Yt| p )t>o for p £ {1, . . . , 7}. 

For j £ N and a £ R+, we denote by W ha the closure of C°°(R,R) with respect to the norm: 

1/2 

\9\\w* :=!£/_ T-T^dx ) , (6.6) 




where g^ is the k th derivative of g. The space W^ ,a endowed with the norm ||.||^> defines a Hilbert 
space. We denote by W~ 3,a the dual space. Let C J ' a be the space of functions g with j continuous 
derivatives and such that 

| 5 ( fc )(x)| 



When endowed with the norm: 



Vfc < 7, lim ; — — = 0. 

~ H-S-+OC 1 + \x\ a 



\9 (k \x)\ 



\\g\\ ci , a :=£sup (6.7) 
k <jx&R l + Fl 

these spaces are Banach spaces, and their dual spaces are denoted by C~ J ' a . 
In the sequel, we will use the following embeddings (see [Tl H5]): 

C-7,0 ^ W 7,l s W 5,2 ^ C 4,2 ^ W 4,3 ^ ^3,3 ^ ^3,4 ^ £,2,4 

C-2,4 ^ w -3,4 ^ C -3,3 ^ W -4,3 ^ £-4,2 ^ ^-5,2 w ~7,l ^ ^-7,0^ (g g) 

where H.S. means that the corresponding embedding is Hilbert-Schmidt (see [T] p. 173). Let us explain 
briefly why we use these embeddings. Following the preliminary estimates of [43J (Proposition 3.4), it 
is possible to choose VF -3,4 as a reference space for our study. We control the norm of the martingale 
part in VF -4 ' 3 using the embeddings T^ 4 ' 3 <— »■ C 3,3 ^ VF 3,4 . We obtain uniform estimate for the norm 
of rjf in C -4,2 . The spaces W~ 5,2 and W~ 7,1 are used to apply the tightness criterion in [35J (see our 
Lemma l6.8p . The space C -7,0 is used for proving uniqueness of the accumulation point of the family 
{r] T )r>o- 
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Proposition 6.4. Let T > 0. The sequence (i] T )t&r + converges in D([0, T],C 7 '°) when T — > +00 
to the unique solution in C([0, T], C~ 7 ' ) of the following evolution equation: 

(mJ)= f / (Lf{x) + Jf{x))r, s {dx)ds + VWWtU), (6.9) 
JO Jr v ' 

where W(/) is a Gaussian martingale independent of W and which bracket is V(f) x t with: 

V(f) =J R ( r J Q 1 (/(?*) + /((! " " f(x)) 2 G(dq) + 2a 2 (x)f(x) 2 ^ n(dx). (6.10) 

Notice that unlike the discrete case treated in [17] . our fluctuation process has here a finite varia- 
tional part. 

6.3 Proofs 

We begin by establishing the evolution equation for r) T that are announced in Proposition 16.11 
Proof of Proposition \6.1\ From Lemma 12.41 and applying (|2.22p with ft(x) = e _rt//2 f(x), we obtain: 
(Z t+T , f) e" r(t+T)/2 =(Z T , f) +M T (f) 

+ fj (Lf{x) + Jf(x)) e~ r < s+T V 2 Z s+T {dx) ds, (6.11) 

where M?(f) is a square integrable martingale with quadratic variation: 

(M T (/)) 4 = / rf* / e~ rs Z s (dx) r / (f(qx) + /((l — g)x) — /(x)] G(dq) + 2a 2 (x)f'{x) 2 , 
Jt Jr l Jo k ' 

which is the bracket announced in (|6.5p . Computing in the same way (Zt, f) e _ri//2 and taking the 
expectation gives, with (|4.2[) and Proposition! 



Q t f(x) e rt / 2 = /(*) + jf* Q s (L/ + Jf) (x) e^ 2 ds. 

Integrating with respect to Zt and multiplying by e~ rT / 2 imply: 

(Z T ,Q t f)e- r ^/ 2 = (Z T ,f)e- rT / 2 + f e+V-V* da (Z T , Q s (Lf + Jf)). (6.12) 

Jo 

We deduce the announced result from (|6.1|) , (|6.11|) and (|6.12|) . ■ 

We now prove that our fluctuation process r] T can be viewed as a process with values in W~ 3,4 , 
by following the preliminary estimates of [33] (Proposition 3.4). This space W~ 3,A is then chosen as 
reference space and in all the spaces appearing in the second line of (|6.8j) that contain W~ 3 ' 4 , the 
norm of rjF is finite and well defined. 

Lemma 6.5. Let T > 0. There exists a finite constant C that does not depend on T nor on T such 
that 

sup Ejlkl^-s.J < Ce r < T+T ) . (6.13) 



te[o,x] 



32 



Proof. Let (</? p ) p sn* be a complete orthonormal basis of W 5, that are C°° with compact support. We 
have by Riesz representation theorem and Parseval's identity: 



^ 9 f ( z t+T,f P ) 2 (Z T ,Q t (p p ) 



E[N 7 



(6.14) 



Under the Assumption 16.21 (iii) and thanks to Remark 14.11 and Example [H we use the same proof as 
in Theorem g21 especially (|4.9|) and (|4.1U|) : 



0<E, 



(Z t+T ,(Pp) 2 (Z T ,Q t ip p ) 2 



E[N t+T } 2 E[N T ] 2 



-r(t+T) 



1 iiQt+TVp + r nQsJ-i [Qt+T-s^p <S> Qt+T-sVp) e" rs ds 
+ e~ rT nQ T {QtWp) 2 + r I nQsh^T-sQuPp® Qr-sQtVp) e_rs ds 

i-t+T rl i- 

<2e~ rT nQt+Tpl + ^r / <pl(qx) e~ rs fiQ t+T (dx) G(dq) ds 

Jo Jo Jr 



(6.15) 



since by (|3.14|) . Cauchy-Schwarz' inequality and symmetry of G: 



J2(Qt+T-sWp\ ® Qt+T-s\<p P \){x) =2 (Q t+ T-s\<Pp\(qx)Qt+T-s\<Pp\{(l - q)x))G(dq) 

<2 J Q t+T _ 8 <p 2 p (qx)G(dq). 

We deduce from flgJH) and (lfU5l) that: 

e-^E^H^H^a] <4e~ rT [ J^*) nQ t+T {d. 

+8 



p> 

t+T r l 



y^ipl(qx)e rs fj,Q t+ T(dx)G(dq)ds. 

10 JO JTt^ 

Let us consider the linear forms D x ^p(g) = g(Fx) for F £ [0, 1], x G R and g £ W 3,4 ^ C 2,4 : 

\D x>F (g)\ = \g(Fx)\ < (1 + |x| 4 )||g|| C 2,4 < C(l + \x\ 4 )\\g\\ w3A 
Using Riesz representation theorem and Parseval's identity, we get: 

D x ,f(Vp) 2 = \\D x ,f\\w-*>* ^ °( 1 + l x l 4 )' 

p>i 



(6.16) 



We deduce from (|4.2|) . (|6.16|) and Assumption 16.21 (iv) that: 



e- rit+T) M\\riJ\\ 2 w-^] <C^ 



l + \Y t 



t+T\ 



1 _ e ~r(t+T) 



<c, 



(6.17) 



(6.18) 



where the constant C is finite and does not depend on T nor T. This completes the proof. 
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We now turn to the proof of the central limit theorem stated in Proposition 16.41 To achieve this 
aim, we first prove the next Lemma on moment estimates. 



Lemma 6.6. Suppose that Assumption I6'.£l is satisfied and let T E M + . 

(i) We have: 

sup supE^ [|| rj[ 1 1 £7-4,2] < +00. 

TeK+ i<T 

(ii) Let us denote by Mj the operator that associates Mj (/) to f. Then 

sup sup E^ [|| Mj Hv^-4,3] < +00. 

TeK+ t<T 



(6.19) 



(6.20) 



Proof. Let us first deal with (|6.2U|) . We consider the following linear forms: D x ^{g) = o~(x)g'(x) and 
D x , q (g) = g(qx) + g((l - q)x) - g(x). Notice that for g £ W 4 < 3 ^ C 3 ' 3 , x GR and q G [0, 1], 



\D x>a (g)\ = \a(x)g'(x)\ < <t(l + \x\)\g'(x)\ < C(l + |x| 4 )|| 5 || C 3, 3 < C(l + |x| 4 )|| ff || W 4, 3 , 
\D*M\ = \9{qx) + - q)x) - g{x)\ < 3(1 + |x| 3 )|| 5 || C 3,3 < C(l + |x| 3 )||g|| W 4,3, 



(6.21) 



where C does not dependent on x nor on q. This implies that D x ^ a and -D^g are continuous from W 4,3 
into M, and their norms in W~ 4,3 are upper bounded by C(l + |x| 4 ) and C(l + |x| 3 ) respectively. Let 
us consider a sequence of functions (</? p ) p eN* constituting a complete orthonormal basis of W 4,3 and 
that are C°° with compact support. Using Riesz representation Theorem and Parseval's identity, we 
get 



p>l 
We have 



lA^H^-4,3 < C(l+|x| 8 ) and J^A^foO 

p>i 



< C(l + |x| 6 ). (6.22) 



E^ [ sup \\Mjf w ^} < E M [ ]T sup M 4 T (v9 p ) 2 



< 



p>l 
r-T+T 



4^E M [(M T (^ p )) T 

Z s (<ix) 



E, 



< C 



c 



T+Y 



T+T 



ds E„ 



E[iV, 

Z s {dx) 
E[N S ) 

U + I^l 8 ) 



(• 



P >i p>i 



(l + |x| 8 ) 



(6.23) 



where the first inequality comes from [T] Lemma 6.52, the second is Doob's inequality, the third line 
is a consequence of (|6.5p . the fourth inequality comes from the bounds (|6.22|) and the last equality 
comes from (|3.1|) . The proof is then finished since by Assumption 16.21 (iv). sup i>0 E^[Y^ 8 ] < 00. 

Let us now consider the proof of (|6.19p . Recall J defined by (|6.3p . It is clear that J is a bounded 
operator from C 4 ' 2 into itself: 

\\J<p\\cv < C\\<p\\cw, (6.24) 
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where C does not depend on ip G C 4,2 . 

Let us denote by C/(i) the semi-group of the diffusion with generator L given by (|5.ip . Proposition 
3.9 in 03] and Assumptions IQ1 yield that for <p G C 4 ' 2 and V G C 3,3 



sup||Z7(t)(^)|| C 4, 2 <C|M| C 4, 2 and sup ||Z7(t)(^)|| C 3, 3 < 



(6.25) 



where C does not depend on ip nor on ip. 

Let us consider the test function ipt '■ ( s > x ) ^ f (i ~~ • s )v 9 ( x ) with </? G C 4 ' 2 . Using Ito's formula: 



{r£,JU(t-s)<p)ds+ / (dM?,U(t-s)(p), 



that is r/f = / J7(i - s)*J*rff ds + [ U(t - s)* dMj ', where U(t - s)* and J* stand for the adjoint 

Jo Jo 
operators of U(t — s) and J. We deduce that for t < T: 



E^llrfH 2 ^] <2T / E^ \\U(t-s)*J* V J\\ 2 c . 4 
Jo 1 

Thanks to (f6T2i|) and (f6T25|) . we have for s < t < T, 



ds + 2E. 



Z7(t -s)*dM s T ||^_ 4 ,2 



E, 



ds < CE,, 



I T 1 1 2 

l 7 7s He- 4 - 2 



ds. 



(6.26) 



(6.27) 



The second term of the r.h.s. of (|6,26[) is upper bounded by considering the norm in W 4 ' 3 . To prove 
that 



sup supE M 

TeR+ t<T 



U(t - s)*dM^ 



< +oo, 



(6.28) 



we use similar arguments as those used for the proof of (|6.2U|) and (|6.25|) . In the proof below, we 
replace the linear forms D Xja and D xs by 4,(_ s , ff and D X)t _ ss with D X)t - Bj(T (<p) = D x ^ a {U(t - s)ip) 
and D Xtt -. St q(<p) = D x , q {U(t - s)ip). Notice that by (RT2TD for g G IU 4 - 3 ^ C 3 ' 3 , i£l and q G [0, 1], 

\D xM {g)\ = \D x>a (U(t)g)\ < C{1 + |x| 4 )||[/(%|| C 3.3 < C(l + |x| 4 )|| 5 || c3 ,3 < C(l + |x| 4 )|M| W 4 >3 , 
lA^G?)! = |i>x )ff (^(^)| < C(l + |x| 3 )||?7(t) 5 || C 3,3 < C(l + Ixl 3 )^!!^, 

where C does not dependent on x. Using again Riesz representation Theorem and Parseval's identity, 
we get 



D x ,tA<Pp) 2 = HAcll^-4,3 < C(l + |x| 8 ) and ]T D xM {ip p f = || A^H^ < C(l + \ x \% 



p>i 



p>i 



(6.29) 
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where C does not dependent on x nor on q. We have with the same arguments as in (|6,23[) : 

E M [sup|| fu{t-s)*dMj\\ 2 w ^] <Ej]Tsup [\u(t- S )<p p ) 2 dM' 
t<Y Jo p> it<rJo 



P>1 

rT+Y 



ds E 



< C 
= C 



T+Y 



ds E 



T+Y 



ds E 




£ Dx,t-,«{<Pp? G(dq) + 2 £ Ar. 
p>i p>i 



(1 + \x\ 



t-s 



<r(<Pp) 



Z s (dx) 



EiV, 



The proof is then done, as sup t>0 E M [Y^ 8 ] < oo by Assumption 16.21 (iv). 
Thus we get from (JEMD, (RT2TD and <KW\ : 



E,[\\r 1 ?\\ 2 c - 4a ]<C[l + 



E, 



o 



I T\\2 

Ws \\c- 4 ' 2 



ds 



We use Gronwall's Lemma and the fact that E^x [| | ??/ H^-4,2] is locally bounded (see Lemma \<i.5\) to 
conclude. ■ 

We now prove the tightness of the fluctuation process. 

Proposition 6.7. Let T > 0. The sequence (?/ t )tgm + is tight in B([0, T], W~ 7,1 ). 

We use a tightness criterion from [38], which we recall (see Lemma C p. 217 in |43j ) . 



Lemma 6.8. (see Lemma C p. 217 in JJ3lj ) 

A sequence (6 T )reM+ °f Hilbert Li-valued cddldg processes is tight in D([0, T],H) if the following 
conditions are satisfied: 

(i) There exists a Hilbert space Hq such that Hq ^h.S. H and \/t < T, sup<p gR+ E[||0^||^ o ] < +00, 

(ii) (Aldous condition) For every e > 0, there exists 5 > and Tq G R + such that for every sequence 
of stopping time tj- < T, 



sup supP(||G 

T>T <;<S 



T 

TT+<, 



e 



\H 



> e) < e. 



Proof of Prop. 16'. 71 We shall use Lemma 16.81 with Hq = W 5,2 and H = W 7,1 . Condition (i) is 



J||2 



< 



a direct consequence of the uniform estimates obtained in (|6.19|) and of the fact that \\rjt \\ w s,2 

r , iin T ii 2 

u W'lt He- 4 - 2- 

Let us now turn to condition (ii). By the Rebolledo criterion (see e.g. |38j). it is sufficient to 
show the Aldous condition for the the finite variation part and for the trace of the martingale part 
of (|6.4p . Let (v?p)p>i be a complete orthonormal system of W 7 ' 1 > C 6,1 . We recall that the trace of 
the martingale part is defined as ti w -7,\((M T )) t = J2 P >i(M' 1 \fp))t (see e.g. Joffe and Metivier [35]). 
Let e > and tt < T be a sequence of stopping times. For Tq > and 5 > 0, following the steps of 



36 



(HEUD, we get: 



sup supP(|tr l4 /-7,i((M T )) TT+? -ti w -7,i((M T )) TT \ > e) 

T>T ?<5 v J 



(6.30) 



< sup sup -E 

t>t ?<5 e 



TT 4[iV s+T ]' r ,/ 



D x>q (tp p ) 2 G{dq) + 2 ^ D x>a (ip p f)ds 
p>i p>i 



Using the embedding W 7 ' 1 > C 6,1 and computations similar to (|6.21|) . we obtain that: 

P>1 

Thus (ET301) gives: 



p>i 



, „,|ii,,- 7 ,i < C(l + |x| 2 ) and D x ^{ip p ) 2 = || Ae^II^-t.i < C(l + |x| 



sup supP M (|tr w/ -7,i((M T )) TT+? -tr iy - 7 ,i((M T )) TT | > e) 
T>T ?<<5 v 7 



<— sup E M 

£ T>T 

C 



\ TTTf T A 7 1 * 1 I 



T T 



s+T 



ds 



4\ „-r(s+r T +T) 



<— sup E M / (Z S+TT+T , 1 + \x\ ) e 
e T>T '-Jo 

<- sup TeJe. [(Z s+r ,l + |x| 4 )e-^ +T ) 

e T>T JO 



ds 
ds 



(6.31) 



by using the strong Markov property of (Zf)t>o- Now, using the branching property: 



E Ztt (Z s+T ,l + |x| 4 )e 



4\ p -r(s+T) 



E,, 



Z TT (dy) 



(Z s+T ,l + \x\*)e- r ^ T \ 
= f E y [l + \Y s+T \ 4 ]Z TT (dy) 

< [ ((tt,1 + \x\ 4 ) + B(l + \y\ 4 )e-^ +T Az TT (dy). 



(6.32) 



for some j3 and B > given by (|5.7|) (see Remark fo. 31 (h)). Since we have a Yule tree, E[iV rT ] < E[A^x] = 
exp(rT). Moreover, using (|2.22|) where the integrand in the second term of the r.h.s. is negative for 
our choice f(x) = |x| 4 and noticing that Z s is a positive measure, we obtain with localizing arguments 
that for any t S R + : 



E M [(Z iA7T ,l + |x| 4 )] <(fi,l + M 4 > + / (86 + 24a)E /1 [(Z aMT ,l + \x\ 4 )}ds. 

Jo 

We deduce from Gronwall's lemma that: 

E^Z^l + M 4 )] <^,l + |x| 4 )e( 8 " b+24 ^ T . 
Then (JOB), (|Q2l) and ([633]) imply that: 



sup supP M (|tr w - 7 ,i((M T }} rr+ ,-tr K ,-7,i((M T )) TT | > e) < — (e rT + e ( 85+24<f)T ) 



T>T ?<<5 



(6.33) 



(6.34) 



which finishes the proof of the Aldous inequality for the trace of the martingale. 
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Remark 6.9. Notice that this also shows that (M t )t>o is tight in W 7>1 . 
For the finite variation part: 

( 

T>T ?<<5 v ' JO JO 



. rr T +c; pt t 

supsupPlll/ (L + J)*r]Jds — / (L + J)*r]Jds\\ w _ 71 > 

T>Tn c<<5 JO JO ' ' 



< sup sup -^-E 

T>T ?<5 £ 

< sup sup -^-E 

T>T ?<5 £ 

C5 rr- 

< sup — 



Iw- 7 ' 1 



(L + J)*?7^||L-r,i^ 



E 



< 



T>T S" JO 

C6(T + 5) 



I T ll 2 

Wt llc-4,2 



eft 



sup supE 

t>t t<r 



I T ll 2 



(6.35) 



We use Cauchy-Schwarz' inequality for the second inequality. The third inequality is obtained by 
noticing that under the Assumption 16.21 and for ip £ W 7 ' 1 : 



\\Ltp\\ C 4,2 < C\\ip\\ C 6,i < C\\(f\\ w r,i 



(6.36) 



as W 7,1 "—^ C ' . We can make the r.h.s. of (|6.35p as small as we wish thanks to (|6.19j) . and this ends 
the proof of the tightness. ■ 

Then, we identify the limit by showing that the limiting values solve an equation for which unique- 
ness holds. This will prove the central limit theorem. 



Proof of Proposition pT^l First of all, by Remark 16.91 the sequence of martingales (M t )t>o is tight 
in W~ 7,1 and thus also in C~ 7 ' by (|6.8p . Let us prove that in the latter space, it is moreover C-tight 
in the sense of Jacod and Shiryaev [37] p. 315. Using the Proposition 3.26 (iii) of this reference, it 
remains to prove the convergence of sup t<T HAM^I 1,^-7,0 to where AM t T = Mf — Mf_. Since the 
finite variation part of (|6.4p is continuous, AM t T = A-qJ , and since in (|6.ip t \-+ (Z t , Qtf) is continuous, 
we have for / G C 7,0 : 



r(t+T) 



sup \AM[(f)\ =supe-^ \f(q(ui,t + T)X^ t+T) ) + /((l - q(u,t + T))X%% 



H+T 



(6.37) 

the label of the particle that undergoes division at t+T, and where q(ui, t+T) 
n the splitting. By convention, if there is no splitting at t + T, the term in 
OTP is 0. Thus su Pi < T |AM t T (/)| < 3e" rT / 2 H/IU < 3e" rT / 2 II flU™. 



where u(u, t+T) £ Vt+T is ^ , , 

is the fraction which appears in the splitting. By conver 
the supremum of the r.h.s. of (|6.37|) is 0. Thus sup t<T 
This proves that: 



C 7,0. 



sup 1 1 AM/ \\c-r,o < 3e 



-rT/2 



(6.38) 



which converges a.s. to when T — > +oo. This finishes the proof of the C-tightness of (M t )t>o in 
C~ 7,0 . The inequality (|6.38p also ensures that the sequence sup t<T ||AM 4 T || 1^-7,1 is uniformly inte- 
grable. From the LLN of Proposition 14.21 the integrand of (|6.5|) converges to W x V(f) which does 
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not depend on s any more. Since W is n e >oO"(^ e )-measurable, it follows that W and W are indepen- 
dent. Thus, using Theorem 3.12 p. 432 in [37], we obtain that (M t )t>o converges in distribution in 
D([0, T], C~ 7,0 ) to a Gaussian process W with the announced quadratic variation. 

By Proposition 16.71 the sequence (?7 T )t>o is tight in W~ 7,1 and hence also in C~ 7,0 by (|6.8p . 
Let tj be an accumulation point in B([0, T], C~ 7 ' ). Because of (|6.4p and ()6.38p . rj is almost surely a 
continuous process. Let us call again by (t] T )t>o, with an abuse of notation, the subsequence that 
converges in law to r/. Since rj is continuous, we get from (|6.4|) that it solves (|6.9|) . Using Gronwall's 
inequality, we obtain that this equation admits in C([0, T], C~ 7 ' ) a unique solution for a given Gaussian 
white noise W which is in C -7 ' . This achieves the proof. ■ 
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